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SECTION  1 
INTRODUCTION 

An  analytical  solution  is  presented  for  the  stresses  and  strains  around  a  circular 
hole  in  a  Mohr-Coulomb  material  under  conditions  of  plane  strain  when  subjected  to 
a'hydrostatic  free-field  pressure.  The  solution  is  a  direct  extension  of  the  solution  of 
Florence  and  Schwer  (Florence,  1978)  and  (Florence,  undated)  to  include  arbi  raiy 
dilatancy  and  the  case  where  out-of-plane  yielding  occurs  in  the  free  fiel  - 

1.1  BACKGROUND. 

The  basic  approach  to  this  problem  was  defined  by  (Newmark,  1969)  in  response  to 
the  need  for  a  method  of  designing  underground  structures  to  resist  the  effects  of 
ground  shock  due  to  nuclear  weapons.  In  the  development  of  his  solution,  Newmar 
assumed  an  elastic-perfecUy-plastic  medium  with  plastic  volume  constancy  m  a 
plane  strain  configuration.  Newmark's  work  was  expanded  by  Hendron  and  Aiyer 
(Hendron,  1972),  who  provided  the  solution  for  a  dilatant  material  using  an  associ¬ 
ated  flow  rule.  Both  the  Newmark  and  the  Hendron  and  Aiyer  solutions  were  based 
on  the  assumption  that  the  out-of-plane  stress  was  the  intermediate  pnncipa 
stress  There  are,  however,  certain  combinations  of  applied  load  and  maten 
properties  where  the  out-of-plane  stress  is  not  the  intermediate  principal  stress. 

Florence  and  Schwer  (Florence,  1978)  and  (Florence,  undated)  generalized  the 
Hendron  and  Aiyer  solution  by  eliminating  the  requirement  that  the  out-of-plane 
stress  be  the  intermediate  principal  stress.  Their  solution  is  considerably  more 

complex,  since  it  requires  one  to  deal  with  multiple  yield  conditions  in  a  sing  e 

material.  In  fact,  the  solution  is  divided  into  two  cases  depending  on  the  material 
properties  and  loads. 

Detoumay  and  others  (Detoumay,  1983)  presented  a  more  general  solution  which 
allowed  arbitrary  dilatancy  through  the  use  of  a  nonassociative  flow  rule,  I bu 
solution  maintained  the  requirement  that  the  out-of-plane  stress  be  the  mteimedi- 
ate  principal  stress.  A  similar  solution  was  developed  independently  by  (Mer  , 
1982),  although  we  have  not  found  it  in  published  form. 
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Wintergerst  and  others  (Wintergerst,  1991)  generalized  both  cases  of  the  Florence 
and  Schwer  solution  to  allow  arbitrary  dilatancy  through  the  use  of  a  nonassociative 
flow  rule,  thus  providing  a  solution  with  both  arbitrary  dilatancy  and  no  limitation 
on  the  intermediacy  of  the  out-of-plane  stress. 

The  analytical  solutions  developed  by  Florence  and  Schwer  and  by  Wintergerst,  et  al 
for  the  plane  strain,  axisymmetric  compression  of  a  circular  opening  in  a  Mohr- 
Coulomb  material  make  the  assumption  that  out-of-plane  yielding  of  the  free -field 
does  not  occur.  These  solutions  assume  that  an  elastic  zone  always  exists  outside 
the  plastic  annulus  and  therefore  do  not  apply  to  the  free-field  yielding  case  (i.e.,  an 
infinite  outer  radius  of  the  plastic  zone).  However,  as  will  be  shown  herein,  certain 
combinations  of  friction  angle,  elastic  properties,  and  free-field  pressure  can  cause 
the  free-field  to  yield  before  significant  deformation  of  the  opening  has  occurred, 
rendering  these  solutions  incomplete.  Yielding  of  the  free-field  does  not  necessarily 
translate  into  excessive  closure  of  the  tunnel. 

1.2  SCOPE  AND  ORGANIZATION  OF  THE  REPORT. 

The  earlier  solutions,  including  the  provision  for  arbitrary  dilatancy,  are  extended  in 
this  report  to  include  yielding  of  the  free  field. 

Section  2  contains  a  complete  problem  description,  a  discussion  of  the  solution 
approach,  the  definitions  of  the  Florence  and  Schwer  solution  cases,  and  a  derivation 
of  the  conditions  under  which  the  free  field  will  yield.  In  developing  the  solution  it  is 
convenient  to  separate  the  problem  into  two  general  cases  based  on  how  the 
concentric  plastic  zone(s)  develops  around  the  opening.  These  two  cases  were 
identified  by  Florence  and  Schwer  as  Case  I  and  Case  II.  The  difference  between 
these  two  cases  is  that  the  initial  plastic  zone  that  develops  at  the  edge  of  the  open¬ 
ing  is  governed  by  two  different  yield  conditions,  depending  on  the  specific 
combinations  of  material  properties  and  loading.  In  Case  I,  the  out-of-plane  stress 
within  the  initial  plastic  zone  is  the  intermediate  principal  stress.  In  Case  II, 
however,  the  out-of-plane  stress  is  the  maximum  (least  compressive)  principal 
stress.  Section  2  also  contains  a  derivation  of  the  material  and  load  parameters 
which  demarcate  Cases  I  and  II.  Since  additional  concentric  plastic  zones  may  form 
outside  of  this  initial  zone,  these  two  cases  are  each  subdivided  into  two  subcases. 
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For  Case  I,  there  are  conditions  that  cause  only  one  initial  plastic  zone  to  form  and 
no  free-field  yielding  is  possible.  These  conditions  constitute  Case  la  solutions. 
Cases  where  additional  plastic  zones  develop  with  continued  loading  will  constitute 
Case  lb  solutions.  The  free  field  may  or  may  not  yield  depending  on  the  materia 
properties.  The  Case  I  solutions  are  described  in  Section  3. 

For  Case  II  the  order  of  the  formation  of  the  plastic  zones  differentiates  the  two 
solution  subcases.  In  Case  Ha,  the  initial  plastic  zone  is  followed  by  the  formation 
of  a  second  plastic  zone,  which  also  begins  at  the  edge  of  the  hole.  Under  certain 
conditions  further  loading  will  cause  the  free  field  to  yield.  In  Case  lib  the  free  field 
yields  after  formation  of  the  initial  plastic  zone.  Further  loading  causes  a  third 
plastic  zone  to  form.  This  third  plastic  zone  has  the  same  stress  conditions  as  the 
second  plastic  zone  in  Case  Ha,  and  it  also  begins  at  the  edge  of  the  hole.  The  Case 
II  solutions  are  described  in  Section  4. 

Tables  summarizing  the  calculational  procedures  applicable  to  the  various  solution 
cases  and  radial  zones  are  presented  in  Section  5.  Numerical  examples  of  the 
extended  solutions  to  the  tunnel  problem  are  presented  in  Section  6.  These  exam¬ 
ples  include  comparisons  with  the  results  of  finite  element  analyses.  Finally, 
references  are  listed  in  Section  7 . 

Frequently  used  symbols  are  defined  in  Table  1-1  which  is  adapted  from 
(Wintergerst,  1991). 
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Table  1-1.  Definitions  of  frequently  used  symbols. 


a 

E 

f 


Pa 

p! 

Pb 

Pb 

Pb 

Pb 


A 

Pb 

A  * 

Pb 

Pb 

r 

R 

R 

R 

R' 

ou 

<*r 

<*0 

er 

£e 

Note: 


v 

♦ 

X 

to 


-  interior  radius  of  the  opening 

-  Young's  modulus  for  medium 

-  yield  surface  definition 

-  plastic  potential  surface  corresponding  to  yield  surface  f 

-  shear  modulus  for  medium 

-  arbitrary  dilatancy  parameter 

-  friction  parameter 

-  internal  pressure  acting  on  hole  boundary 

-  negative  of  the  radial  stress  at  the  elastic-plastic  boundary 

-  far-field  pressure  (compressive  stress)  at  large  radius 

-  far-field  pressure  at  initial  yield 

-  far-field  pressure  at  which  R  =  R 

-  far-field  pressure  when  the  inner  plastic  zone  begins  to  form  under 

Case  II  yield  conditions 

-  far-field  yield  pressure 

-  pressure  value  used  to  discriminate  Cases  Ila  and  lib 

-  far-field  pressure  at  which  the  third  plastic  zone  forms  in  Case  lib 

-  arbitrary  radius  to  any  point  in  the  medium 

-  radius  to  the  elastic-plastic  boundary 

-  maximum  radius  at  which  the  radial  and  out-of-plane  stresses  are  equal 

-  minimum  radius  at  which  the  radial  and  out-of-plane  stresses  are  equal 

-  radius  to  the  elastic-plastic  boundary  when  pb  =  pb 

-  unconfined  compressive  strength  of  medium 

-  stress  in  radial  direction 

-  stress  in  circumferential  direction 

-  stress  in  out-of-plane  direction 

-  strain  in  radial  direction 

-  strain  in  circumferential  direction 

-  strain  in  out-of-plane  direction 

Strains  are  further  distinguished  by  the  superscripts  (e)  to  denote 
elastic  component  and  (p)  to  denote  plastic  component,  and  overdots 
to  denote  derivatives  with  respect  to  time 

-  Poisson's  ratio  for  medium  (0  <  v  <  0.5) 

-  friction  angle  for  medium 

-  flow  constant  of  proportionality 

-  dilatancy  angle 
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SECTION  2 

PROBLEM  DESCRIPTION  AND  APPROACH 

The  purpose  of  this  section  is  to  furnish  a  complete  description  of  the problem  and 
the  approach  taken  to  solve  it.  In  addition,  we  present  derivations  wluch_(l) 
the  two  solution  cases  developed  by  Florence  and  Schwer;  and  (2)  define  the 
renditions  under  which  the  free  field  will  yield.  These  derivations  are  i^udedin 
this  section  since  they  are  a  part  of  the  ffee-field  yielding  solutions  for  both  Case 
and  Case  II  presented  in  Sections  3  and  4,  respectively.  The  problem  geomet^ 
loading,  material  law,  and  solution  restrictions  are  described  in  Section  2.1.  The 
overall  solution  approach  is  discussed  in  Section  2.2.  The  Florence  and  Schwer  Case 

I  and  Case  II  conditions  are  developed  in  Section  2.3.  Finally,  t  econ  ions  or 

free-field  yielding  are  derived  in  Section  2.4. 


2.1  PROBLEM  DESCRIPTION. 


The  problem  geometry  consists  of  a  circular  hole  with  radius  a  in  an  infinite  medium 
under  conditions  of  plane  strain.  The  loading  of  the  material  consists  of:  (1)  a 
compressive  pressure,  Pl,  applied  at  the  inside  edge  of  the  hole  (r  equals  a),  and  (2) 
a  compressive  pressure,  pb,  applied  at  r  equal  to  infinity  (define  as  e  ee 
Tensile  values  of  p.  and  p„  are  precluded.  The  loading  is  assumed  to  be  monotomc, 
with  p.  equal  to  pb  until  p.  reaches  its  final  value.  Subsequently,  p.  is  held 
constant  while  pb  is  increased  to  its  final  value.  The  internal  pressure,  p  can  be 

considered  to  represent  the  resistance  provided  by  a  structural  finer.  The  na  v  ue 

of  p  would  then  represent  yielding  of  the  liner.  The  loading  is  ass  umed  to  be 
aoulied  sufficiently  slowly  that  inertial  effects  can  be  ignored.  Finally,  it  is  assumed 
that  during  the  initial  loading  phase  with  p.  equal  to  p„,  the  materia  response  1 

linear. 

The  material  is  assumed  to  be  governed  by  a  Mohr-Coulomb  elastic-plastic  consti- 
tutive  law  with  a  failure  surface  defined  by 
f  =  a3-Na1  +  ou 

where  o,  is  the  maximum  (least  compressive)  principal  stress  and  c3  is  the 


5 


and  pb,  are  positive  in  compression.  The  material  parameter  au  represents  the 
unconfined  compressive  strength  of  the  material.  The  material  parameter  N  is 
related  to  the  friction  angle,  <>,  by 

N  =  l+sin(<|>)  (2.2) 

l-sin(<t>) 

The  previous  solutions  developed  for  this  problem  are  not  valid  for  N  equal  to  one, 
that  is,  it  has  been  assumed  that  the  friction  angle  was  greater  than  zero.  We  make 
the  same  assumption  for  the  solution  presented  herein. 


Plastic  flow  of  the  material  is  assumed  to  be  governed  by  the  nonassociative  flow 
rule 


(2.3) 


where  represents  the  components  of  the  plastic  strain  rate  tensor,  A ,  is  the 
plastic  multiplier  (a  scalar),  g  represents  the  plastic  potential  surface,  and  CJ^ 


represents  the  components  of  the  stress  tensor.  The  form  of  the  plastic  potential 
surface  is  taken  to  be 


g  =  a3-Mcr1  +  K 

where  the  material  parameter  M  is  related  to  the  dilatation  angle  co  by 

l+sin(co)  (2.5) 

l-sin(a>) 

and  k  is  an  arbitrary  constant.  The  constant  k  is  arbitrary  since  only  the  deriva¬ 
tives  of  g  with  respect  to  the  stresses  are  used  in  the  flow  rule.  Note  that  if  co  is 
equal  to  <j> ,  then  M  is  equal  to  N  and  the  flow  rule  becomes  associative. 

Furthermore,  the  value  of  M  is  restricted  by 


1<  M  <  N 


(2.6a) 


or 


0  <  co  <  <)) 

A  dilatancy  angle  of  zero  (M  =  1)  indicates  a  non-dilatant  material. 


(2.6b) 


It  should  be  noted  that  for  the  geometry  and  loading  of  this  problem,  the  principal 
stresses  coincide  with  the  normal  stresses  or,  ae,  and  cz,  where  the  subscripts  r,  6, 
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and  z  represent  the  radial,  circumferential,  and  out-of-plane  components, 
respectively.  The  shear  stress  components,  ort ,  a0z,  and  are  a  zero. 

Elastic  material  response  is  assumed  to  be  isotropic  and  therefore  governed  by  two 
independent  elastic  constants.  Throughout  this  report  we  use  Poisson  s  ratio,  v. 
Young's  modulus,  E,  and  shear  modulus  G,  depending  on  which  is  most  convenient. 
The  relationship  between  G,  E,  and  v  is  given  by 

E  (2-7) 

G  =  2(I+v) 

We  further  restrict  the  value  of  Poisson's  ratio  to  be  greater  than  zero  and  less  than 
one-half. 


2.2  APPROACH. 


The  general  solution  strategy  used  by  Florence  and  Schwer  for  nonassociative  flow 
(Florence,  1978)  and  subsequently  adopted  by  (Wintergerst,  1991)  for  arbitrary 
dilatancy  is  as  follows.  The  medium  is  broken  up  into  regions  (zones)  which  depend 
on  the  stress  state  (both  principal  stress  ordering  and  magnitudes)  and  hence  on  the 
governing  yield  condition  (if  any).  The  number  of  regions,  and  the  stress  state  m 
each  region  is  dependent  on  both  the  loading  and  material  properties.  The  stress 
states  include  both  those  governed  by  elastic  material  response  and  those  goveme 
by  plastic  material  response. 


Florence  and  Schwer  used  numerical  analysis  to  guide  their  postulates  about  the 
existence  and  conditions  of  the  various  plastic  zones.  After  deriving  them  analytical 
solutions,  they  performed  finite  element  analyses  and  obtained  excellent  agreement 
between  analytical  and  numerical  results.  We  have  adopted  the  same  solution 
philosophy  in  extending  the  analytical  solution  to  include  yielding  of  the  free  field. 


Within  each  region,  equilibrium  conditions,  strain-displacement  compatibility,  an 
the  constitutive  law  (including  yield  conditions  and  flow  rules,  where  applicable) 
were  used  to  derive  an  ordinary  differential  equation  (in  terms  of  a  primary 
unknown)  which  governs  the  response  in  that  region.  Depending  on  the  region  the 
primary  solution  variable  was  taken  to  be  either  the  radial  displacement  or  e 
radial  stress.  The  choice  was  based  on  ease  of  formulation  and  solution.  The  other 
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unknowns  (i.e.,  the  remaining  stress  and  strain  components,  or  the  displacement) 
are  written  in  terms  of  the  primary  unknown. 

A  general  solution  for  the  governing  differential  equation  in  each  region  is  then 
determined.  These  solutions  typically  contain  unknown  constants,  which  are 
determined  either  by  applying  continuity  conditions  at  the  region  boundaries  or  by 
using  applied  loading  conditions  at  the  edge  of  the  hole  or  at  r  equal  to  infinity. 

2.3  CASE  I  AND  CASE  II  DEFINITIONS. 

In  this  section,  the  Florence  and  Schwer  derivations  of  the  Case  I  and  Case  II 
conditions  (Florence,  1978)  and  (Florence,  undated)  are  described. 

The  derivation  of  the  elastic  solution  for  a  hole  in  an  infinite  medium  under  condi¬ 
tions  of  plane  strain  is  well  documented,  so  only  the  results  are  presented  here. 


Or  =-Pb+(Pb-Pa)^2- 

(2.8) 

<*e=-Pb-(Pb-Pa)^2- 

(2.9) 

= -2vpb 

(2.10) 

a2 

2Ger  =-(l-2v)pb+(pb-pa)^j 

(2.11) 

a2 

2Ge0  =  — (1— 2v)pb  —  (pb  Pa )  ~~2 

(2.12) 

For  pa  equal  to  pb  during  the  initial  loading,  Equations  2.8  and  2.9  reduce  to 

af  =  ae  _  _Pb.  For  values  of  Poisson's  ratio  satisfying  the  inequality  0  <  v  <  0.5,  the 

stress  ordering  is  or  =  ae  <  oz.  With  increased  loading,  yielding  of  the  material  will 


be  governed  by 

f1  =  oe-Naz+au  =0 
f2  =  Cr-NCTz+<V=0 

Substitution  of  the  stress  components  in 
Equation  2.13a  gives 


(2.13a) 

(2.13b) 


terms  of  pb  into  the  yield  condition  of 
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(2.14) 


_  _  u 

Pa  "  Pb  “  1_2NV 

For  compressive  values  of  p„  (equal  to  p.),  the  material  will  not  yield  if  2Nv  >  1. 

For  the  condition  2Nv  <  1.  the  material  will  yield  throughout  when  Equation  2.14  is 
satisfied  As  discussed  in  Section  2.1,  we  assume  that  the  internal  pressure,  p„ 
does  not  cause  yielding  of  the  material,  and  thus  we  restrict  the  internal  pressure  to 

c. 


Pa  < 


when  2Nv  <  1 


(2.15) 


l-2Nv 

So  far,  there  is  no  restriction  on  p,  for  2Nv  >  1,  but  one  trill  be  developed  later. 


As  discussed  in  Section  2.1,  once  the  internal  pressure  reaches  its  maximum  value, 
it  is  held  constant,  while  the  free-field  pressure  is  increased.  Once  the  external 
pressure  is  increased  above  the  internal  pressure,  we  have  (according  to  Equataons 
2.8  through  2.10)  «,  <  o,  and  <J„  <  <T,  throughout  the  medium.  At  the  edge  of  the 

hole  (r  equals  a),  ot(a)>o,(a)  if  p.  <  2vpk  and  o,(a)  <  o.(a)  if  P.  >  2vpb.  An 
inspection  of  Equations  2.8  through  2.10  shows  that  the  maximum  stress  difference 
occurs  at  the  edge  of  the  hole,  and  therefore  initial  yielding  will  occur  there. 


We  will  denote  the  external  pressure  at  which  yielding  begins  as  p„.  The  stresses  at 
the  edge  of  the  hole  (r  equals  a)  when  yielding  commences  are  evaluated  using 


Equations  2.8  through  2.10. 

CTr(a)  =  “Pa 

<*e(a)  =  -2Pb+Pa 

az(a)  =  _2VPb 


(2.16a) 

(2.16b) 

(2.16c) 


An  inspection  of  Equations  2.16  indicates  that  if  p.  <  2vp„,  then  the  stress  order  is 

0  <  a  <  cr  at  r  equals  a,  and  therefore  the  governing  yield  condition  is 

"  A  (2.11 

f  =  ae-Nar  +  au=0 

Substitution  of  the  stresses  from  Equations  2.16  into  the  yield  condition  of 
Equation  2.17  leads  to  the  initial  yield  pressure 


Pb=f[(N+l)Pa+°u] 

Substitution  of  Equation  2.18  for  pb  into  the  inequality  pa  <  2vpb  gives 


(2.18) 
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(2.19) 


Pa  < 


VG„ 


when  (N  +  l)v  <  1 


1-(N  +  l)v 

The  restriction  (N+  l)v  <  1  is  imposed  since  we  are  only  considering  compressive 
values  of  pa.  For  (N  +  l)v  >  1  there  is  no  restriction  on  pa.  The  condition 
ae  <  az  <  ar  at  initial  yielding  constitutes  Case  I. 

Inspection  of  Equations  2.16  also  indicates  that  if  2vpb  <Pa  <Pb>  then  the  stress 
order  at  the  edge  of  the  hole  is  ce  <  or  <  Gz  and  therefore  the  governing  yield  condi¬ 
tion  is  given  by 


f  =  oa-Nc  +o„  =0 


(2.20) 


Substitution  of  Equations  2.16  into  Equation  2.20  leads  to  the  initial  yield  pressure 

p.  =  Pa+G-^-  when  Nv  <  1  (2-21) 

Pb  2(1- Nv) 

The  restriction  Nv  <  1  is  imposed  since  we  are  only  considering  compressive  values 
of  pb.  The  condition  ae  <  aT  <  cz  at  initial  yielding  constitutes  Case  II. 

For  the  condition  Nv  >  1,  yielding  governed  by  Equation  2.20  will  not  occur  for 
compressive  values  of  Pb.  Instead,  as  pb  is  increased  it  will  eventually  reach  a  value 
where  pa  <  2vpb  so  that  the  stress  ordering  is  oe  <  oz  <  ar  and  yielding  is  governed 
by  Equation  2.17.  Thus,  for  Nv  >  1,  we  have  a  Case  I  condition  and  therefore  pb  is 
given  by  Equation  2.18.  Also,  for  the  condition  Nv  >  1,  it  follows  that  (N  +  l)v  >  1  so 
no  restriction  is  placed  on  pa  (as  was  done  in  Equation  2.19). 

For  Case  II  with  Nv  <  1,  the  yield  pressure  of  Equation  2.21  is  substituted  into  the 
inequality  2vpb  <  pa  <  pb  (for  which  ce  <  ar  <  oz  at  the  edge  of  the  hole)  which  gives 


Pa+^L. 


2v(Pa  +  <V)  < 

2(1-Nv)  Pa  2(1-Nv) 


Equation  2.22  can  be  rearranged  to  give 


vo„ 


<  Pa  < 


- - —  when  0  <  Nv  <  — 

l-2Nv  2 


l-(N+l)v 

- - <  p  when  —  <  Nv  <  (N  +  l)v  <  1 

1-(N+  l)v  Pa  2 

Equation  2.23b  provides  the  restriction  on  pa  for  2Nv  >  1  alluded  to  in  the 
discussion  of  Equation  2.15. 


(2.22) 

(2.23a) 

(2.23b) 
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In  summary,  the  conditions  under  which  initial  yielding  is  governed  by  the  yield 
condition  Equation  2.17  constitute  Case  I.  The  condit.ons  under  whrch  mitia 
yielding  is  governed  by  the  yield  condition  of  Equation  2.20  constitute  Case  II. 
the  physical  difference  between  Case  I  and  Case  II  is  the  stress  ordering  at  the  edge 
of  the  hole  when  yielding  begins.  The  Case  I  and  Case  II  conditions  are  summarized 
in  Table  2-1  (adapted  from  Table  I  of  (Florence,  1978)). 

2.4  FREE-FIELD  YIELDING  CONDITIONS. 


(2.24) 


Equilibrium  in  the  free  field  dictates  that 

i  e.  the  in-plane  principal  stresses  are  equal  to  the  negative  of  the  externally 
applied  pressure.  The  negative  sign  comes  from  the  convention  that  stresses  axe 
taken  positive  in  tension,  while  the  applied  pressures  (pa  and  Pb)  are  taken  posit 
in  compression.  For  elastic  conditions,  the  out-of-plane  stress  oz  may  be  obtained 

from  Hooke's  law,  the  plane  strain  condition,  and  Equation  2.24. 

,  I.  n  (2.25a) 

Eez  =  cz  -  v(or  +  oe)  =  0 

(2.23b) 

cz  = ~2vpb 

where  E  is  Young's  modulus  and  e,  is  the  out-of-plane  strain  (zero  by  definition). 


For  compressible  materials  (v  <  0.5),  the  order  of  the  principal  stresses  is 


c  =  cfl  <  Cz 


(2.26) 


(2.27a) 

(2.27b) 


which  leads  to  failure  (yielding)  being  governed  by  the  intersection  of  the  following 
two  failure  surfaces. 

f1  =  or-Ncz  +  cu=0 
f  2  =  ce  -  Ncz  +  cu  =  0 

The  magnitude  of  the  applied  pressure  at  initial  yield  in  the  free  field  (pb )  is 
determined  by  substituting  Equations  2.24  and  2.26  into  either  of  Equations  2.27a 

or  2.27b, 

-pb+2Nvpb+cu  =0 
or 


(2.28a) 
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Table  2-1.  Case  I  and  Case  II  conditions. 


Property  Relations 

Internal  Pressure 

Case 

Nv  <  (N  +  l)v  <  1 

-  va« 

Pa  1-(N+  l)v 

I 

1  <  Nv  <  (N  +  l)v 

or 

Nv  <  1  <  (N  +  l)v 

Pa  >0 

I 

—  <  Nv  <  (N  +  l)v  <  1 

2 

V<Tu 

Pa  1-(N+  l)v 

II 

0  <  Nv  <  - 
2 

van  .  <JU 

l-(N+l)v  a  "  l-2Nv 

II 
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(2.28b) 


l-2Nv 


Note  that  if  1-  2Nv  <  0,  then  free-field  yielding  can  only  occur  for  tensile  values  of 
p“scussed  in  Section  2.1,  tensile  values  of  pb  and  p.  are  precluded. 
Therefore,  for  l-2Nv  <  0,  the  earlier  solution  (Wintergerst,  1991)  is  complete. 
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SECTION  3 
CASE  I  SOLUTION 

In  this  section,  the  arbitrary  dilatancy  solution  for  Case  I  is  extended  to  include  the 
effects  of  free-field  yielding.  It  was  shown  in  Section  2.2  that  yielding  of  the  free 

field  occurs  when 

Pl>  g»—  =  p.  for  l-2Nv > 0 
Pb-l-2Nv  F 

Although  the  earlier  solution  (Wintergerst,  1991)  is  valid  for  Pb  <  Pb,  the  entire  Case 
I  solution  is  presented  below  for  completeness. 

As  stated  in  Section  2.1,  pb  is  assumed  equal  to  pa  (the  internally  applied  pressure) 
until  pa  reaches  its  final  value.  The  external  pressure,  pb,  is  then  increased 
monotonically  while  pa  is  held  constant.  The  initial  material  response  is  assumed 

to  be  elastic  (Figure  3-1).  As  discussed  in  Section  2.3,  initial  yielding  of  the  material 
begins  at  the  edge  of  the  hole  (r  equals  a)  when  pb  =  [(N  +  l)pa  +  cu]/2,  with  stress 
conditions  ce  <  crz  <  oT,  and  is  governed  by  the  failure  surface 

f  =  g0  -  Nor  +  au  =  0  (3'2) 

After  initial  yielding,  an  inner  plastic  zone  forms  with  outer  boundary  radius  R. 
Outside  of  R,  the  material  is  elastic  (as  shown  in  Figure  3-2).  For  certain  combi¬ 
nations  of  material  properties,  az  is  always  the  intermediate  principal  stress 
(throughout  the  medium  and  regardless  of  the  value  of  pb);  thus,  no  additional 
plastic  zones  will  form  and  the  free  field  will  not  yield.  In  this  report,  the  conditions 
under  which  az  is  always  the  intermediate  principal  stress  will  be  termed  Case  la. 
For  the  case  where  the  plastic  flow  rule  is  associative  (M  equal  to  N),  the  Case  la 
solution  is  the  Hendron  and  Aiyer  solution  (Hendron,  1972). 

The  Case  la  solution  can  be  considered  to  consist  of  two  phases. 


Phase  1.  Elastic  solution 

Phase  2.  One  plastic  zone  with  an  outer  elastic  zone 

For  other  combinations  of  material  properties,  oz  is  not  always  the  intermediate 
principal  stress  throughout  the  medium.  As  pb  increases,  additional  plastic  zones 
can  form  where  oz  is  the  maximum  principal  stress.  As  pb  increases  further,  the 
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free  field  may  yield  if  l-2Nv>0.  The  conditions  where  0,  is  not  always  the 
intermediate  principal  stress  constitute  Case  lb.  The  demarcation  between  Case  la 
and  Case  lb  conditions  is  discussed  in  Section  3.1. 


For  Case  lb,  the  initial  inner  plastic  zone  increases  in  size  as  loading  continues, 
until  its  boundary  reaches  a  limiting  value  R,  where  o,(R)  =  o,(R).  This  inner 
plastic  zone  boundary,  defined  by  r  =  R ,  remains  fixed  as  p.  is  increased.  Continue 
loading  causes  the  formation  of  two  more  plastic  zones.  In  the  middle  plastic  zone, 
the  stress  condition  is  o9  <  a.  =  ur  and  yielding  is  governed  by 


f1-o,-Not  +  <J.-0 
f2  =  o,-Na,  +  o„  =  0 

In  the  outer  plastic  zone,  the  stress  condition  is  <  or  <  az 


(3.3a) 

(3.3b) 

and  yielding  is  governed 


by  the  single  failure  surface 

(3.4) 

f  =  g6-Nc2  +  gu  =  0  _ 

The  boundary  between  the  middle  and  outer  plastic  zones,  defined  by  r  -  R,  where 
a  (R)  =  oz(R).  Outside  of  the  plastic  zones,  the  material  response  was  assumed  in 

the  earlier  solutions  to  be  elastic.  The  boundary  between  the  outer  plastic  zone  and 
the  elastic  zone  is  defined  by  r  =  R.  The  three  plastic  zones  and  the  elastic  zone, 
with  their  boundaries  and  stress  conditions  are  shown  in  Figure  3-3. 


With  further  loading,  the  free  field  may  yield  if  1-  2Nv  >  0.  If  the  free  field  yields,  a 
fourth  plastic  zone  (Figure  3-4)  is  postulated  to  exist.  In  subsequent  discussions, 
the  fourth  zone  is  designated  the  "far  outer  plastic  zone".  This  term  is  used 
throughout  this  report  to  indicate  free-field  yielding,  even  though  there  may  be  as 
few  as  two  plastic  zones.  The  stress  condition  in  this  zone  is  aT  -  oe  <  a2  an  e 


yield  conditions  are 

f  x  =  ar  -  Naz  +  ou  =  0 
f2  =  ae— Naz+au  =° 


(3.5a) 

(3.5b) 


Based  on  the  progression  of  the  formation  of  plastic  zones,  the  solution  for  Case  lb 
can  be  considered  to  consist  of  four  phases: 

Phase  1.  Elastic  solution 

Phase  2.  One  plastic  zone  with  an  outer  elastic  zone 
Phase  3.  Three  plastic  zones  with  an  outer  elastic  zone 
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Figure  3-3.  Phase  3  conditions  for  Case  lb 


Phase  4.  Four  plastic  zones  (free  field  yielding) 

Note  that  Phase  1  and  Phase  2  are  identical  for  Cases  la  and  lb. 

The  elastic  solution  was  presented  in  Section  2.3.  The  solutions  for  the  phases  with 
plastic  zones  are  presented  in  Sections  3.1  through  3.3. 

3.1  CASE  la  AND  lb,  PHASE  2  SOLUTION. 

In  this  section,  we  develop  the  solution  for  Phase  2  of  Case  I.  As  discussed  previ¬ 
ously,  in  Phase  2  there  is  a  single  (inner)  plastic  zone  with  an  outer  elastic  zone 
(Figure  3-2).  The  derivation  presented  below  closely  follows  those  contained  in 
(Florence,  1978)  and  (Wintergerst,  1991).  The  demarcation  between  Cases  la  and  lb 
is  also  defined. 

Initial  yielding  of  the  material  occurs  at  the  edge  of  the  hole  (r  =  a)  when  pb  is  equal 
to  pb  (derived  below).  As  discussed  previously,  for  Case  I  conditions  Ge  <  Gz  <  gt  and 

the  yield  condition  is 

f  =  CJe  -  Norr  +  <JU  =  0  (3'2) 

Substitution  of  the  elastic  solution  for  ar  and  ae  (presented  in  Section  2.3)  evalu¬ 
ated  at  r  equals  a  into  the  yield  condition  of  Equation  3.2  leads  to  the  following 
expression  for  the  initial  yield  pressure,  pb . 

pl  =  i[(N+1)p.+<t.]  (3'6) 


The  equation  of  equilibrium  is  given  by 

d<rr 


r— r-^  +  Or  -C7e  =  0 


(3.7) 


dr 


and  is  valid  throughout  the  material.  Using  the  yield  condition  of  Equation  3.2,  ce 
can  be  written  in  terms  of  ar, 

a,=N<J,-<Ju  (3'8) 

so  the  equilibrium  Equation  3.7  can  be  rewritten  as 

r^-<rr(N-l)+au=0  (3*9) 

dr  , 

The  solution  to  Equation  3.9,  with  the  boundary  condition  cr( a)  =  -pa,  is  given  by 
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Pa  + 


N-lJU 


N-l 


+ 


Ou 

N-l 


(3.10) 


The  strain  rates  are  determined  by  the  nonassociative  flow  rule,  discussed  in 
Section  2.1.  For  the  yield  condition  of  Equation  3.2,  the  plastic  strain  rate  compo- 

nents  are  given  by 


e(p)  =  X—  =  -XM 
r  dcr 


dca 


=  x^-=o 


(3.11a) 


(3.11b) 


(3.11c) 


da 


(3.13a) 

(3.13b) 


where 

(3.12) 

g  =  ce  -  Mar  +  k 

and  X  is  the  plastic  multiplier.  Since  the  loading  is  assumed  to  be  monotonic  and 
the  material  law  is  rate  independent,  the  overdot  in  Equations  3.11  may  be  in  er- 
preted  to  represent  either  differentiation  with  respect  to  time  or  differentiation  wi 
respect  to  pb.  From  the  flow  rule  of  Equations  3.11  it  follows  that 

e(p)+Me(ep)  =  0 

e<zp>=0 

which,  when  integrated  with  respect  to  time,  give 

.  .  ,  ,  N  (3.14a) 

e(rp)+Me'p)  =h1(r) 

,  .  ,  ,  ,  (3.14b) 

e(zp  =h2(r) 

The  strain  expressions  in  Equations  3.14,  evaluated  at  a  fixed  radius  r,  remain 
constant  as  loading  proceeds  (i.e.,  they  are  independent  of  time,  or  eqmvalently  of 
load  Pb).  When  the  elastic-plastic  radius  R  is  initially  at  this  radius  r,  the  plastic 

strains  are  zero,  which  implies  that  hx(r)  =  0  and  h2(r)  -  0,  therefore 

,  ,  (3.15a) 

e(p)  +  Me(ep)  =  0 

(3.15b) 

elp)=0 

The  strain  rates  are  assumed  to  be  decomposed  into  elastic  and  plastic  parts 


21 


er  =  i[e>  +  e(rP> 


P  -  P(e)  +  £(p) 

te  —  t-e  T  fce 


(3.16a) 
(3.16b) 

ez=E(ze)  +  e(zp)  (3-16c) 

Integrating  Equation  3.16a  with  respect  to  time  gives  £r  -(e(re)  +  e(rP>)  =  h(r)-  When  r 

equals  R,  the  strain  is  entirely  elastic  so  er  =  e(re)  and  e(rp)  =  0  giving  h(r)  =  0.  Using 
the  same  argument  for  the  other  two  strain  components  leads  to 

(3.17a) 
(3.17b) 
(3.17c) 


er  =  e(re)  +  e(rp) 


£0  =  ele)  +  £ 


(p) 


ez  =  e(ze)  +  e'z 

The  plane  strain  condition  requires  that  ez  =  0,  and  the  flow  rule  led  to  £zp  -  0, 
therefore  £(,e)  =  0 . 


,(p) 


The  elastic  strains  are  related  to  the  stresses  by  Hooke's  law. 


Ee(re)  =ar-v(ae  +  cz) 

(3.18a) 

E£(6e)  =a9-v(ar  +  az) 

(3.18b) 

Eeze)  =  az  -v(ar  +  a„) 

(3.18c) 

From  Equation  3.18c  and  the  condition  £ze)  =  0 

cz=v(cr+ae) 

(3.19) 

Substituting  Equation  3.19  for  az,  and  Equation  3.8  for  a6  into  Equations  3.18a 
and  3.18b  leads  to 

2Ge(re)  =[l-(N+l)v]ar+vau 

2Ge(0e)  =  [N  -  (N  +  l)v]ar  -  (1-  v)ou 

where  ct  is  given  by  Equation  3.10. 


(3.20a) 

(3.20b) 


The  strain-displacement  relations  (for  small  strains)  are  given  by 

du 
dr 


er 


e9=- 


u 

r 


(3.21a) 

(3.21b) 


where  u  represents  the  radial  displacement.  The  compatibility  equation  is  derived 
by  eliminating  u  from  Equations  3.21. 


22 


3L+E,-er  =  0 

dr 


(3.22) 


Using  the  strain  decomposition  of  Equation  3.17,  Equation  3.15a  to  eliminate  e',"1, 
Equations  3.20  for  the  elastic  strain  components,  and  Equation  3.10  for  a„  t  e 
compatibility  Equation  3.22  can  be  rewritten  as 

r^+(N+l)er  =^(N  +  l)[(N-l)P.  +  ^](f)  <3.23) 

The  solution  to  Equation  3.23,  with  the  boundary  condition  e(ep)(R)  =  0  is 

.I  1-v  N  +  l  r/xr  (3.24) 


=  _  I"  V  N  +—  [(N  -  l)pa  +  CTU  ]  - 
2G  M  +  Nlv  '  JU 


The  remaining  unknown  is  the  elastic-plastic  boundary  radius  R.  Using  R  and  p  in 
phice  of  a  and  p.,  respectively,  in  the  elastic  solution  given  in  Section  2.3  (Equations 

2.8  through  2.10)  gives 

(3.25a) 

ar(R)  =  -p 

(3.25b) 


oe(R)  =  -2pb  +  p 


(3.25c) 


az(R)  =  -2vpb 

This  notation  is  used  throughout  the  report.  That  is,  when  stress  or  strain  com- 
ponents  are  immediately  followed  by  a  single  variable  in  parentheses,  that  vanab  e 
denotes  the  spatial  location  where  the  stress  or  strain  component  is  being 
evaluated.  For  example,  the  notation  or(R)  denotes  the  stress  function  cr  evalu¬ 
ated  at  r  equal  to  R.  This  rule  will  also  apply  to  the  radial  displacement  function,  u. 


The  yield  condition  of  Equation  3.2  still  holds  at  r  equals  R,  so 

lr,  x  l  (3-26) 

Pb=-|[(N  +  1)P  +  a»] 

G  (TU-  2Pfc~q-»-— p  (3'27' 

°r(R)  N+1  P 

Setting  r  equal  to  R  in  Equation  3.10  and  equating  the  result  to  Equation  3.27  gives 
f  a  ¥Rf  ~X  ■  -  2pb~gu  (3.28a 

ip-  +  N^lJuJ  N-l  N  +  l 


ou  2pb-au 


(3.26) 


(3.27) 


(3.28a) 


N-l 


N  +  l 
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BY-1  =  2  (N-l)pb  +  au  (3.28b) 

a )  N+l  (N-l)pa+au 


The  solution  presented  above  is  valid  with  continued  loading  as  long  as  oz  <  or  in 
the  plastic  zone.  From  Equations  3.19  and  3.2 

or-oz  =  [l-(N+l)v]or  +  vou  (3- 

In  Equation  3.29,  oz  <  or  when  the  right  hand  side  of  the  equation, 

[l-(N  +  l)v]or  +  vou,  is  greater  than  zero.  If  (N+l)v>  l,then  [l— (N  +  l)v]  is 
negative.  Since  or  is  compressive  (negative)  and  vou  is  positive,  then 
[l— (N  +  l)v]or  +  vou  is  positive  and  thus  oz  <  or.  If  (N  +  l)v  <  1,  then 

[  1  -  (N  + 1)  v]or  +  vou  is  positive  for  or  <  -  +  ^  summary '» 


(1)  oz  <  or  when  (N  +  l)v  >  1 

vo 

(2)  oz  <  or  when  (N  +  l)v  <  1  and  or  <  ~1_^N"  ^ 

When  (N  +  l)v  >  1,  the  solution  presented  above  is  complete  since  no  more  plastic 
zones  will  form  and  the  free  field  will  not  yield.  This  case  will  be  referred  to  as  Case 
la.  In  the  second  case,  when  (N  +  l)v  <  1,  the  minimum  value  of  or  occurs  at  r  equal 
to  R,  so  that  as  loading  proceeds  the  elastic-plastic  boundary  radius  R  attains  a 
limiting  value,  denoted  R ,  where 


vo., 


and 


a.oa-a.ao—i^ 


(l-v)g. 


(3.30) 


(3.31) 


(N+l)v 


The  magnitude  of  the  loading  when  R  equals  R  is  denoted  as  pb  which  is  given  by 
Equations  3.25,  3.26,  and  3.30 

-  q,  (3.32) 

Pb  2[l-(N  +  l)v] 

and  the  radius  of  the  elastic-plastic  interface  can  be  determined  from  Equations 
3.281)  and  3.32. 


va  V 


(l-2v)ou 


[l  -  (N  +  l)v][(N  -  l)Pa  +  ou  ] 


(3.33) 
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The  relationship  for  R  is  independent  of  pb,  as 
will  remain  constant  with  increased  loading. 


is  the  stress  field,  and  therefore  both 


3.2  CASE  lb,  PHASE  3  SOLUTION. 

The  solution  for  Phase  3  of  Case  lb  is  developed  in  this  section.  As  discussed  previ¬ 
ously,  there  are  three  plastic  zones  and  an  outer  elastic  zone  in  Phase  3  (Figure  3-3). 
Again,  the  derivation  presented  below  closely  follows  those  contained  in  (Florence, 

1978)  and  (Wintergerst,  1991). 

Florence  and  Schwer  postulated  that  with  increased  loading  (pb  >Pb)>tw0  addi' 
tional  plastic  zones  are  formed  outside  of  the  initial,  inner  plastic  zone  In  the  outer 
plastic  zone  (R  <  r  <  R)  the  stress  condition  is  o„  <  <5,  <  a„  therefore  the  yield 

condition  is 

„  (3.4! 

f  =  ce-Na2  +  ou=0 


The  plastic  strain  rates  are  obtained  from  the  nonassociative  flow  rule. 
e(rp)  =  0 
e(9p)  =  ^ 


(3.34a) 

(3.34b) 


e'p)  =  -MX, 

Using  similar  arguments  for  integrating  the  plastic  strain  rates  as 
Section  3.1  leads  to 


(3.34c) 
presented  in 


e(p)  +  M4p)=0  (3.35a) 

e(P)  =  0  (335b) 

Using  the  strain  decomposition  of  Equation  3.17,  the  plane  strain  condition,  and 
Equations  3.35  leads  to  the  following  strain  expressions. 


e.  =  e(re) 


(3.36a) 


£e  ~  ee  +  M  z 


e  =0 


(3.36b) 

(3.36c) 
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Using  Hooke's  law  (Equations  3.18),  the  strain  Equations  3.36,  and  the  yield  condi¬ 
tion  of  Equation  3.4  gives  the  following  expressions  for  the  radial  and  tangential 
stress  components  in  terms  of  the  radial  and  tangential  strain  components. 

[MN  + 1-  (M  +  1)(N  +  l)v]cr  =  [MN  + 1-  (M+ N)v]2Ger 

+M(N+l)v2Ge6+(M-1)vau  (3.37a) 

[MN  + 1-  (M  +  1)(N  +  l)v]ae  =N(M  +  l)v2Ger 

+MN2Gee-[l-(M  +  l)v]au  (3.37b) 


Substitution  of  the  stress  Equations  3.37  into  the  equilibrium  Equation  3.7,  and 
using  the  strain-displacement  Equations  3.21,  leads  to  the  following  ordinary 
differential  equation  in  terms  of  the  radial  displacement  u, 

q2u__1_(±^Kt  (3.38) 

P  MN  2G 

where 


d2u  du 


dr2 


■  +  r 


dr 


(M-N)v  2 
MN  H 


MN 


R2  = - 

p  MN+1-(M  +  N)v 

The  solution  to  Equation  3.38  is 
2Gu  =  A^1  +  A2r~Y2  +  Ar 

where 


=_iizMa=-d- 


l-2Nv0u 


(l-2v)pb 


Yi 


Y2 


(N-M)v  ,  RN^M)V 

P  on/rvr  P  1  /oiv/rNTl2 


2MN 


(2MN)  P2 


(M-N)v  2 


2MN 


i 


(M-N)  v2  |  1 
(2MN)2  +P2 


(3.39) 


(3.40) 


(3.41) 


(3.42) 


(3.43) 


and  Ax  and  A2  are  constants  that  are  yet  to  be  determined.  In  Equation  3.41,  the 
definition  of  pb  from  Equation  2.28b  is  used.  The  strains  can  be  computed  from  the 

strain-displacement  relations  and  the  displacement  Equation  3.40. 

2Ger  =  YiA^1-1  -  Y2A2r“Y2_1  +  A  (3.44a) 

2Gee  =  A^1'1  +  A2r_T2_1  +  A  (3.44b) 

The  strain  Equations  3.44  can  be  substituted  into  the  stress  Equations  3.37  to  yield 
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Or=CrlA1rT‘"1+Cr2A2r'72_1-Pb 
ae  =  CelA1rY,'1  +  C02A2r"T2"1-Pb 

where 

Crl  =  -{[MN  + 1-  (M  +  N)v]y!  +  M(N  +  l)v} 

c 

Cr2  =  — {-[MN  +  1-(M+N )  v]y2  +  M(N  + 1)  v} 

c 

Cei  =  ^[N(M  +  l)vYa  +  MN] 

Ce2  =  i[-N(M+  1)vy2  +MN] 

C  =  MN  + 1-  (M  +  1)(N  +  l)v 


(3.45a) 

(3.45b) 

(3.46a) 

(3.46b) 

(3.46c) 

(3.46d) 

(3.46e) 


At  this  point  the  unknowns  are  the  two  constants  A1  and  A2  and  the  plastic  zone 
boundary  radii  E  and  R.  Two  of  the  conditions  used  to  evaluate  the  unknowns  are 
provided  by  enforcing  continuity  of  radial  stress  and  radial  displacement  (or, 
equivalently,  circumferential  strain)  at  the  boundary  with  the  outer  elastic  zone  (r 
equals  R).  The  elastic  zone  stresses  evaluated  at  r  equals  R  can  be  substituted  into 

the  yield  condition,  of  Equation  3.4  to  give 

cr  (R)  =  -2(1-  Nv)pb  +  ou  (3‘47 

The  circumferential  strain  at  r  equals  R  is  given  by 

2G€e  (R)  =  -2(N  -  l)vpb  -  au  v3'48 


Continuity  of  radial  stress  at  r  equals  R  leads  to 
CrlA1RTl_1  +  Cr2A2R_Y2_1  =  2(l-Nv)APb 
while  continuity  of  circumferential  strain  at  r  equals  R  leads  to 
A^1-1  +  A2R'Y2_1  =  2(N  -  l)vAPb 


(3.49a) 

(3.49b) 


where 


(3.50) 


Equations  3.49  constitute  a  system  of  two  equations  in  the  two  unknowns  AjRTl  1 
and  A2R_Yz_1j 


Crl  Cr2"lf  A.R71'1 1  ={2(l-Nv)APbl  (3.51) 

1  1  J{a2R_Y2-1J  \2(N-l)vApbJ 
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which  can  be  inverted  to  give 

JBJ  f  Aft*'1  ]  =  1 

|bJ  [a2e-,,'1J  c„-ct! 

where  Bx  and  B2  are  constants. 


-Cr2"|f2(l-Nv)Apbl 
Crl  J[2(N-l)vApb  J 


(3.52) 


Florence  and  Schwer  postulated  that  at  the  inner  radius  R  of  the  outer  plastic  zone, 
cr(R)  =  cz(R),  so  that  the  yield  condition  of  Equation  3.4  may  be  written  as 

Nor(R)-oe(E)  =  ou  _  (3'53 

Using  the  stress  expressions  in  Equation  3.45  evaluated  at  R  and  substituted  into 
Equation  3.53  leads  to 

(NC„  -  C»1)A1R'"-1  +  (NC„  -  CM  )A2R-’'-'  -  N(l-  2v)p„  =  0  (3.54 

Substituting  the  results  of  Equation  3.52  into  Equation  3.54  leads  to  the  following 
nonlinear  equation 


VRy 


rVY2_1 


(NCrl  — Cei)B1  £  +  (NCr2-Ce2)B2  £  -N(l-2v)pb  =  0 


(3.55) 


which  can  be  solved  for  the  radius  ratio  (R/R) .  Determination  of  the  individual 
values  of  R  and  R  requires  a  description  of  the  state  in  the  middle  plastic  zone. 


As  discussed  earlier,  Florence  and  Schwer  postulated  the  existence  of  a  middle 
plastic  zone  |R  <  r  <  R  j  where  oe  <  crr  =  (Jz,  so  that  yielding  is  governed  by 


fl  =  <70-N<Tr+au  =0 

f  2  =  °e  ~  N°z  +  °u  =  0 

Solving  for  ae  in  terms  of  cr  using  Equation  3.3a,  the  equilibrium  Equation 


(3.3a) 

(3.3b) 
3.7  has 


the  solution 


N-l 


1-1  i 


N-l 


(3.56) 


N-l 


Substituting  R  and  or( R)  from  Equations  3.30  and  3.33,  respectively,  into 
Equation  3.56  gives 


or  =  -  Pa  + 


N-lJVa 


-  I  + 


N-l 


(3.57) 
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Note  that  Equation  3.57  is  identical  to  Equation  3.10,  so  the  radial  stress  is  repre 
sented  by  the  same  formula  throughout  the  combined  inner  and  middle  plastic 
zones  a  <  r  <  R .  Providing  continuity  of  radial  stress  at  the  interface  r  equals  R 


leads  to 

CrlBx 


v  R  y 


+  Cr2B2 


^rVY2_1 

VRy 


-  _  L  ,  gu 

Pb"  [Pa  +  N_! 


vay 


N- 1 


(3.58a) 


or 


_y.-i  N(l-2v)pb-(N-l)C„B,^|j  +(N-l)Ct2B2[R 


TJ  A  Yj'‘ 


R 

a 


_  -  (3.58b) 

(N-l)p.+o„ 

Using  Equation  3.58b,  R  can  be  determined  since  we  have  previously  solved  for  the 
ratio  R/R  from  Equation  3.55  and  the  constants  B,  and  B2  from  Equation  3.52. 
After  solving  for  R,  R  can  be  computed  from  the  ratio  R/R  and  Ax  and  A2  can  be 

computed  using  Equation  3.52. 


The  plastic  strain  rates  are 
e(p>  _  1  dgi  ^  3g2_ 

£'  ~K3g,  !3o. 


p(P> 


~(p) 


,x  Oh. 


doe 


n  dgl  .5,  ^2. 

=  ^1  — +  ^2  d 


da 


given  by  the  nonassociative  flow  rule. 

=  -MX,! 

=  ■+■  X-2 

=  -MX,2 


so  that 

[e(rp)  +  MEqp)  +  e(zp)  ]  =  0 

which,  when  integrated  with  respect  to  time,  gives 
[e(p)  +  Me(ep)  +  e(2p)  ]  =  h(r) 


(3.59a) 

(3.59b) 

(3.59c) 

(3.60) 

(3.61) 


We  designate  the  plastic  strains  at  r  equals  R  as  e,p),  £eP),  and  ez  ,  so  that 

[e>r'>+Mer  +  ^p,]  =  h(K)  _  <3'62) 

But  in  the  outer  plastic  zone  (and  therefore  at  r  equals  R),  the  plastic  strain  rela¬ 
tions  of  Equations  3.35  hold  so  h(R)  must  be  zero.  Since  every  point  in  the  middle 

plastic  zone  will  at  some  time  be  at  R,  then  h(R)  is  equal  to  zero.  Therefore, 
Equation  3.61  becomes 
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(3.63) 


[e(rp)  +  MEeP)  +  e(zp)  ]  =  0 
The  plane  strain  condition,  £z  =  £ze)  +  £(zp)  =  0  >  combined  with  Equation  3.63  gives 

er  =  -Mee  +  Me(0e)  +  £(re)  +  £(ze)  ^ 

With  cr  equal  to  az,  and  oe  eliminated  via  the  yield  condition  of  Equation  3.3a, 

Hooke's  law  can  be  reduced  to 


Ee(re)  =  E£ze)  =  [l  -  (N  +  l)v]or  +  vau  (3-65 

Ee0e)  =  (N  -  2v)ar  -  au  (3'66 

where  cr  is  given  by  Equation  3.57.  Substitution  of  the  elastic  strains  from 
Equations  3.65  and  3.66  into  the  radial  strain  expression  of  Equation  3.64,  and 
using  the  resulting  expression  for  Er  in  the  compatibility  Equation  3.22  leads  to  the 


following  differential  equation. 


2G— (rM*'£  ^_MN  +  2-2(M  +  N+1)v[(N-l)p,+gu] 
dr1  e*J  (N-l)(l+v)  l'  ’  1 


„M+N-1 


.N-l 


(M  +  2)(l-2v)ou  (3.67) 

(N-l)(l+v) 

This  differential  equation  is  directly  integrable.  Integration  over  the  limits  from  r  to 
R  results  in  the  following  expression  for  £e. 


2G£e=2G£e(R) 


^  R  ^M+1 


v  r  y 


,  MN  +  2-2(M+N  +  l)v r/XT _ lb 
(M  +  N)(N-l)(l+v)  LV 


au] 


^r^m+n 


(M  +  2)(l-2v)qu 
(N  -  1)(M  + 1)(1+  v) 


The  circumferential  strain  £e(R)  is  obtained  from  Equation  3.44b. 


(3.68) 


The  radial 


strain  is  obtained  from  Equations  3.64  through  3.66  giving 


2G£r  =  -2GMee  + 


MN  +  2-2(M  +  N  +  l)vg 
1+v 


M-2v 

1+v 


(3.69) 


For  applied  pressures  greater  than  pb,  the  stresses  in  the  inner  plastic  zone  are  still 
given  by  Equation  3.10  for  Gr,  Equation  3.8  for  Ge,  and  Equation  3.19  for  Gz.  The 
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differential  Equation  3.23  for  the  plastic  strain  £?>  still  holds  with  increased 
loading,  but  now  must  be  integrated  from  r  to  R ,  instead  of  from  r  to  R,  whtch 
results  in 


*r\(*r 


2Ge«p)  -  2Ge,p>(ft)l  — 


M  +  N  Va 


(3.70) 


The  other  plastic  strain  components  are  still  given  by  Equations  3.15.  In 

3.70,  e<pl  (r)  is  found  by  setting  r  equal  to  ft  in  Equation  3.68  to  obtain  e, (R) ;  then 

setting  r  equal  to  ft  in  Equation  3.66,  and  using  Equation  3.30  for  ot  (r)  ,  to  obtain 

tf  (ft);  and  then  using  the  strain  decomposition  given  in  Equation  3.16b. 

The  solution  developed  above  breaks  down  when  the  external  pressure,  pb,  reaches 
the  free-field  yield  pressure,  pb.  We  demonstrate  this  as  follows.  From  Equation 
3.50,  as  pb  approaches  pb,  Apb  approaches  zero.  ^ 

lim  Apb  =  lim  (pb  -pb) =  0 

Pb^Pb  Pb->Pb 

Substituting  Equation  3.71  into  Equation  3.52  leads  to 

fO  '70o'I 


lim  AjR*1”1  =  0 

Pb^Pb 

lim  A2R_Til‘1  =  0 

Pb  “*Pb 


(3.72a) 

(3.72b) 


Pb  “*Pb 

It  can  be  shown  that  both  Yl  and  Y2  are  greater  than  aero.  From  Equation  3.72b,  as 
Pb  approaches  pb,  either  A2  goes  to  zero  or  R  goes  to  infinity,  or  both.  Physically, 
one  can  argue  that  as  the  free-field  pressure  approaches  p„,  the  outer  elastic  zone 

vanishes,  with  the  result  that  R  must  go  to  infinity. 

As  pb  approaches  p„,  the  circumferential  strain  at  r  equals  R  is  computed  from 
Equation  3.48  as 


lim  2Gee  (R)  =  -2(N  -  l)vpb  -  au  =  A 


(3.73) 


Pb“>Pb 


where  A  =  -(1_  2v)pb  as  defined  in  Equation  (3.41). 
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From  the  strain-displacement  relations,  it  follows  that 
lim  2Gu(R)  =  AR 

Pb-»Pb 

Equation  3.40  for  radial  displacement  in  the  outer  plastic  zone,  evaluated  at  R, 


gives 

2Gu(R)  =  AiR’'+A2R-’-+AR  <3'75) 

From  Equations  3.74  and  3.75  it  follows  that  as  pb  approaches  pk>  then 

lim  [  AkRYl  +  A2R"Tt  ]  =  0  (3'76) 

Using  Equation  3.72b,  we  showed  that  either  A2  goes  to  zero  or  R  goes  to  infinity,  or 
both.  This  implies  that  Equation  3.76  may  be  decomposed  as  follows 


lim[A2R-H  =  0 

Pb“»PbL  J 

lim  [ A1RYl  1  =  0 

Pb^Pb1  J 

Equation  3.77b  implies  that  A1  must  be  zero  since  Yi  is  greater  than  zero. 


(3.77a) 

(3.77b) 


For  pb  equal  to  pb  (and  Ax  equals  zero),  Equation  3.45  can  be  rewritten  as 

^r=Cr2A2r'Y2"1-Pb 

ae  =  Ce2A2r-^-x-pb 

Again,  we  postulate  that  at  the  inner  radius  R  of  the  outer  plastic  zone, 
ar(R)  =  o,(R),  so  Equation  3.53  still  holds.  Using  the  stress  expressions  in 

Equation  3.78  evaluated  at  R  and  substituted  into  Equation  3.53  leads  to 

(NCr2-C02)A2R"Y2_1-(N-l)pb  =  ou 


(3.78a) 

(3.78b) 


(3.79) 


which  can  be  solved  for 


*2  2  NCr2  -Ce2 


(3.80) 


Continuity  of  radial  stress  at  the  boundary  with  the  middle  plastic  zone  (r  equals 
R)  using  Equations  3.57,  3.78a  and  3.80  gives 


Cr2B2-pb 


r  +  o 

f'R'l 

N-l 

+ 

lP*  +  N-lJ 

UJ 

N-l 

(3.81a) 


or 


va/ 


N(l-2v)pb  -(N-l)Cr2B2 
Pa(N-l)  +  au  “ 


(3.81b) 
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After  solving  for  R,  A2  can  be  recovered  from  Equation  3.80.  The  remainder  of  the 
solution  developed  for  p„  <  p„  still  applies  for  p„  equal  to  pb. 

3.3  CASE  lb,  PHASE  4  SOLUTION. 

The  solution  for  Phase  4.  which  contains  free-field  yielding,  is  presented  in  this 
section.  Note  that  for  Nv  >  I.  the  free  field  will  not  yield  and  therefore  tins  phase  of 

the  solution  is  not  applicable. 

After  the  free  field  yields  (Pb  >  p„),  we  postulate  that  a  far  outer  plastic  zone 
(R  <  r  <  °°)  exists  with  conditions 

V  '  (3.82) 

as  show  in  Figure  3-4.  This  assumption  was  based  on  the  results  of  numerical 
studies  using  this  same  Mohr-Coulomb  material  model.  The  stress  condition 
far  outer  plastic  zone  is  o,=c9<o„so  the  yield  condition  is  given  by 

(3.5a) 

f1  =  ar-Na2  +  ou  =0 

(3.5b) 

f2  =  ce-Ncz+cu=0 

The  expressions  for  stress  in  the  outer  plastic  zone,  discussed  in  Section  (3.3)  are 
repeated  below 

.  .  -y..— i  a  (3.45a) 

cr  =  CrlAxrir  +  Cr2A2r  -  pb 

i  .  _Y2_x  .  (3.45b) 

<ye  =  CelA1rl  +Ce2A2r  Pb 

At  this  point  our  unknowns  are  At,  A2,  R, and  E.  As  discussed  earlier,  E  is  the 
radius  to  the  boundary  between  the  middle  and  outer  plastic  zones.  At  r  equal  to  E , 
it  is  reasonable  to  postulate  that  oe(E)  =  -p„  in  the  outer  plastic  zone.  This 
condition,  along  with  continuity  of  <T,  at  R ,  results  in  two  equations  for  two 
unknowns,  AjR71  1  and  A2R  T2  » 

A1RYl_1l=  |Apbl  <3-83) 

2r_Y2_1J  lApbJ 

which  can  be  inverted  to  give 


Crl  ^r2 
LC01  Ce2j 
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iB  1  f  A  RYl_1 

|Bj  =  A,1R--M"CrlCe2-Cr2C61 


Ce2 

-c 


-c 


ei 


r2 


' rl 


,Apb 
l  Apt 


(3.84) 


where  Ba  and  B2  are  constants.  After  solving  for  B1  and  B2,  the  stress  expressions 

of  Equation  3.45  can  be  written  as 

/  _  yi-1  / 


+  Cr2B2 


o8  =  C81B,5 


NYi- 


+  C02®2 


r  ^  72  1  - 

ij  -Pb 

r  ^  Y2  1  - 

CrJ  -Pb 


(3.85a) 

(3.85b) 


As  discussed  in  Section  3.2,  at  the  boundary  between  the  middle  plastic  zone  and 
outer  plastic  zone  we  have 

Nc,(k)-o8(R)  =  c„  (353) 

Substitution  of  the  stress  expressions  from  Equations  3.85  into  Equation  3.53 
results  in  the  following  nonlinear  equation 


Rvrl 


(NC^-CJB,!  j  I  +(NCr2-Ce2)B2 


'rVY2_1 

R 


-Npb(l-2v)  =  0 


(3.86) 


which  can  be  solved  for  the  radius  ratio  R/R . 


As  developed  in  Section  3.3,  the  radial  stress  in  the  middle  plastic  zone  is 


o¥  r^ 

|_Pa  +  N-  ljva  J 


V-1  a 


N- 1 


Continuity  of  radial  stress  at  the  boundary  R  provides 


CrPl 


1TMT 


Pb  _Pa  +  N-l. 


^  R  ^N_1 


v  a  y 


N-l 


(3.57) 


(3.87a) 


or 


R 


V 


_  N_i  N(l-2v)pb-(N-l)CrlB1  - (N -  l)Cr2B2 

”  vjw _ 


R 

a 


^rVY2_1 

a, 

R 


(N-l)pa  +ou 


(3.87b) 


from  which  we  can  solve  for  R,  since  B1}  B2,  and  R/R  have  been  determined 
previously.  Once  B1(  B2,  and  R  are  determined,  the  solution  presented  in  Section 
3.2  is  valid  for  the  inner,  middle,  and  outer  plastic  zones. 
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In  the  far  outer  plastic  zone,  Hooke's  law  can  be  rewritten  using  Equation  3.82  to 
obtain 

Ee(re)=-pb(l-v)-v<yz 

Ee(ee)=-pb(l-v)-voz 

Ee(zc)  =  az-  2vpb 

The  plastic  strain  rates  are  obtained  from  the  nonassociative  flow  rule. 


(3.88a) 

(3.88b) 

(3.88c) 


e'">  =  X  I|Si  +  k2|S2.  =  -M(X1  +  k!) 

z  da,  da. 


(3.89a) 

(3.89b) 

(3.89c) 


(3.90) 


(3.91) 


which  implies  that 

M(e(rp)  +  e(ep))  +  elP)=0 

which,  when  integrated  with  respect  to  time,  gives 

M(e(rp)  +  £(p) )  +  £(ZP>  =  h(r) 

However,  the  plastic  strain  relations  of  Equation  3.35  still  hold  in  the  outer  plastic 
zone;  hence  at  r  equals  R,  h(r)  must  be  zero.  Therefore,  Equation  3.91  can  be 

rewritten  as 

£<">  =  -Mfe?’  +  4P>)  <3'92) 


From  the  yield  condition  of  Equation  3.5,  the  out-of-plane  stress  is  given  by 


o,=^(-pb  +  <0 


(3.93) 


From  the  plane  strain  condition,  strain  decomposition,  Equations  3.92,  3.93  and 
3.88c 

Eef  =  -Ee‘p>  =  EM(e?’  +  e?1 )  =  -  ^  ^ + 2vPt  <3'94) 

Therefore,  the  sum  of  the  radial  and  tangential  plastic  strains,  £?>  +e‘/>,  is  a  con¬ 
stant  throughout  the  far  outer  plastic  zone,  i.e.,  there  is  no  spatial  variation  in  that 
sum.  We  will  denote  that  constant  as  C. 
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(3.95) 


e(rp)+4p)  EM|  N+N+2VPbI  2GMN(1+ v) 


_g„-(l-2Nv)pb 


The  compatibility  Equation  3.22,  using  strain  decomposition,  can  be  written  as 


,  de<p) 

dr  dr 


+  e(ee)  +  e‘p)-e(re)-£(rP)=0 


(3.96) 


j  (0) 

From  Equations  3.88a  and  3.93,  ^2-  is  equal  to  zero  ,  and  e“  is  equal  to  so  the 


dr 

compatibility  equation  reduces  to 

r^+e(p>_e(rp)  =  0 

dr 

Using  Equation  3.95  to  eliminate  e(rp)  leads  to 
de(p) 


dr 


— +  2eip  -C  =  0 


which  can  be  solved  to  give 


eip)  =^  +  Aa 


--2 


(3.97) 


(3.98) 


(3.99) 


A 

where  A3  is  a  constant.  The  postulation  that  ae(R)  =  -pb  in  the  outer  plastic  zone 

and  the  yield  conditions  of  Equations  3.4  and  3.5  lead  to  all  three  normal-stress 
components  being  continuous  across  the  plastic-plastic  boundary  at  r  =  R .  From 
Hooke's  law,  this  implies  that  ej,e>  is  also  continuous  at  r  =  R.  Continuity  of 
tangential  strain  at  r  =  R  then  implies  that  e(ep)  must  also  be  continuous  at  r  =  R . 
We  invoke  this  condition  to  determine  the  constant  A3 . 


4pI(r)  =  §+a3K-2 


(3.100a) 


or 


A,  =  RN  ee 


>(r)_!|  (3.100b) 

where  e(ep)  (r)  is  computed  from  the  solution  in  the  outer  plastic  zone  (Section  3.2). 
Equation  3.44,  evaluated  at  R,  gives  the  total  strain  £6(r).  Equations  3.45  and  3.4 
can  be  used  to  compute  crr  (r)  ,  ce  (r)  and  o2  (r)  .  After  computing  the  stresses, 
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Hooke's  law  can  be  used  to  compute  the  elastic  strain  4  )(®)- 
strain  e(ep)(R)  is  computed  by  subtracting  e(ee)(R)  from  ee(R). 


Finally,  the  plastic 
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SECTION  4 
CASE  II  SOLUTION 


In  this  section  the  arbitrary  dilatancy  solution  for  Case  II  is  extended  to  include  the 
effects  of  ffee-field  yielding.  Just  as  in  the  Case  I  solution,  the  earlier  solution 
(Wintergerst,  1991)  is  valid  only  if  the  material  response  in  the  free  field  is  elastic. 
Again,  yielding  of  the  free-field  occurs  when 

Pl>_Hs» —  =  p.  when  (l-2Nv)>0  t4-1 

Pb"l-2Nv  Pb  V 

As  with  Case  I,  the  earlier  solution  is  valid  for  pb  <  pb  but  in  the  interest  of  com¬ 
pleteness,  the  entire  solution  is  presented  herein. 


As  derived  in  Section  2.3,  initial  yielding  for  Case  II  begins  at  the  edge  of  the  hole  (r 
equals  a)  with  stress  condition  c6  <  cr  <  o2  and  is  governed  by  the  failure  surface 


f  =  Ga  -  No,  +  ou  =  0 


(4.2) 


This  initial  plastic  zone  increases  in  size  as  loading  continues  until  ar  is  equal  to  az 
at  the  edge  of  the  hole.  At  this  point  a  second  plastic  zone  begins  to  form.  This 
second  plastic  zone  also  begins  at  the  edge  of  the  hole,  and  like  the  first  plastic  zone, 
increases  in  size  as  loading  proceeds.  The  initial  (now  outer)  plastic  zone  is  still 
governed  by  Equation  4.2.  The  second  (inner)  plastic  zone  is  governed  by  the 
following  yield  conditions 

f,  =  ce-Ncr +  cu  =  0  (4'3a) 

f2  =  a,-No,+<5u=0  _  (4'3b) 

The  boundary  between  the  inner  and  outer  plastic  zones,  denoted  as  R,  is  defined  by 


ar(R)  =  o2(R)  ^ 

Outside  of  the  plastic  zones,  the  material  response  was  assumed  in  earlier  solutions 
to  be  elastic.  The  earlier  solutions  can  therefore  be  considered  to  consist  of  three 
phases  (in  order  of  occurrence): 


Phase  1.  Elastic  solution. 

Phase  2.  Elastic-plastic  solution  with  one  plastic  zone  and  an  outer  elastic 
zone. 

Phase  3.  Elastic-plastic  solution  with  two  plastic  zones  and  an  outer  elastic 
zone. 
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In  the  extended  solution,  free-field  yielding  can  occur  either  before  or  after  the 
formation  of  the  second  plastic  zone.  However,  it  is  assumed  never  to  occur  Prior  to 
the  formation  of  the  first  plastic  zone.  This  assumption  is  consistent  with  the 
assumption  that  the  application  of  the  internal  pressure,  p.  (with  p.  equal  to  p„  > 
does  not  cause  material  yielding.  The  case  where  free-field  yielding  occurs  after  e 
formation  of  the  second  plastic  zone  will  be  denoted  as  Case  Ila,  while  the  case 
where  free-field  yielding  occurs  before  the  formation  of  the  second  plastic  zone  wi 
be  denoted  as  Case  nb.  The  demarcation  between  Cases  Ila  and  lib  will  e 
discussed  in  detail  in  Section  4.1. 

When  the  free  field  yields,  a  far  outer  plastic  zone  is  postulated  to  exist.  The  stress 

condition  in  this  zone  is  ar  =  ce  <  oz  and  the  yield  conditions  are 

A  (4.5a) 

f,=ar-Noz  +  au  =0 

A  (4.5b) 

f2  =  ae  -Ncz  +  au  =  0 

This  outer  plastic  zone  occupies  the  domain  R  <  r  < -.with  the  boundaiy  between 
the  initial  plastic  zone  and  the  far  outer  plastic  zone  given  by 

/M  /a\  (4.6) 

^r(R)  =  ae(R) 


Based  on  the  progression  of  formation  of  the  plastic  zones,  the  solution  for  Case  Ila 
can  be  considered  to  consist  of  the  following  four  phases: 

Phase  1.  Elastic  solution. 

aase2.  Elastic-plastic  solution  with  one  plastic  zone  and  an  outer  elastic 


zone. 


Phase 


3.  Elastic-plastic  solution  with  two  plastic  zones  and  an  outer  elastic 


zone. 


Phase  4.  Plastic  solution  with  three  plastic  zones.  The  elastic  zone  is 

eliminated  due  to  free-field  yielding  (i.e.,  it  is  replaced  by  a  far  outer 

plastic  zone). 

For  Case  lib,  the  solution  process  consists  of  the  following  four  phases: 

Phase  1.  Elastic  solution.  . . 

phage  2.  Elastic-plastic  solution  with  one  plastic  zone  and  an  outer  elastic 

zone. 
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Phase  3.  Plastic  solution  with  two  plastic  zones.  The  outer  elastic  zone  of 
Phase  2  is  eliminated  due  to  free-field  yielding  (i.e.,  it  is  replaced 
by  a  far  outer  plastic  zone). 

Phase  4.  Plastic  solution  with  three  plastic  zones. 

Phases  1,  2,  and  4  are  each  identical  for  Cases  Ila  and  lib;  only  the  conditions  in 
Phase  3  differ.  The  conditions  for  each  of  the  phases  are  depicted  in  Figures  4-1 
through  4-5. 

The  elastic  solution  (Phase  1  for  both  Cases  Ila  and  lib)  was  presented  in  Section 
2.3.  The  Phase  2  solution  is  identical  for  Cases  Ila  and  lib,  and  is  presented  in 
Section  4.1.  The  Phase  3  solution  for  Case  Ila  is  presented  in  Section  4.2.  The 
Phase  3  solution  for  Case  lib  is  presented  in  Section  4.3.  Finally,  the  Phase  4  solu¬ 
tion  (identical  for  Cases  Ila  and  lib)  is  presented  in  Section  4.4. 

4.1  CASE  Ila  AND  lib,  PHASE  2  SOLUTION. 

In  this  section,  the  solution  is  presented  for  Phase  2  (initial  yielding  with  one  plastic 
zone  and  an  outer  elastic  zone  as  shown  in  Figure  4-2).  Following  that,  the 
conditions  are  derived  for  the  formation  of  the  second  plastic  zone.  As  will  be  shown, 
these  conditions  also  define  Cases  Ila  and  lib.  Finally,  the  solution  for  when  pb  is 
equal  to  pb  for  Case  lib  is  presented.  The  solution  presented  closely  follows  both 
those  of  (Wintergerst,  1991)  and  (Florence,  undated),  but  with  slightly  different 
notation. 

Initial  yielding  of  the  material  occurs  at  the  edge  of  the  hole  when  pb  is  equal  to  pb 
(determined  below).  As  discussed  earlier,  for  Case  II  conditions,  ce  <  or  <  cz  at  the 

edge  of  the  hole,  and  therefore  the  yield  condition  is 

f  =  o„-N<s,  +  <j„=0  (4'2) 

Substitution  of  the  elastic  solution  for  cr  and  az  (presented  in  Section  2.3)  evalu¬ 
ated  at  r  =  a  into  the  yield  condition  of  Equation  4.2  leads  to  the  following 
expression  for  the  initial  yield  pressure,  pb : 

-  P„+gu  (4.7) 

Pb  2(1- Nv) 
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The  plastic  strain  rates  are  given  by  the  nonassociative  flow  rule  (discussed  in 
Section  2.1) 


where 


e(p>  =  X^-  =  0 

(4.8a) 

r  ter 

k[p)  =  X—  =  X 

(4.8b) 

tee 

e(p)  =  X  dg  =  -MX 

(4.8c) 

tez 

g  =  ae-Maz  +  K 

(4.8d) 

and  X  is  the  plastic  multiplier.  As  discussed  in  Section  3.1,  the  overdots  in 
Equations  4.8  may  be  interpreted  as  either  differentiation  with  respect  to  time  or 
differentiation  with  respect  to  pb  since  monotonic  loading  is  assumed.  From  the 

flow  rule  of  Equations  4.8  it  follows  that 

(4.9a) 


e'p)+Me'p)  =0 
e(p)  =  0 


(4.9b) 


which,  when  integrated  with  respect  to  time,  give 

e<p)+Me<p)  =h1(r)  (4.10a) 

e<p)=h2(r)  (4'10b) 

The  strain  expressions  in  Equations  4.10  evaluated  at  a  fixed  radius  r  remain  con¬ 
stant  as  loading  proceeds  (i.e.,  they  are  independent  of  load  pb).  When  the  elastic- 

plastic  radius,  R,  is  initially  at  this  radius  r,  the  plastic  strains  are  zero,  which 
implies  that  hx(r)  =  0  and  h2(r)  =  0,  and  Equations  4.10  reduce  to 

e(zp)+Me<p)  =  0 

e(p)  =  0 


(4.11a) 

(4.11b) 


The  strain  rates  are  decomposed  into  elastic  and  plastic  parts. 
er  =  e(re)  +  e(rP) 


ee  =  e(ee)+4p) 

ez  =  £(ze)  +  e(2p) 


(4.12a) 

(4.12b) 

(4.12c) 
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Integrating  Equation  4.12a  with  respect  to  time  (or  pb)  gives  er  (e(re)  +  erP  )  h(r)> 

but  when  r  equals  R,  the  strain  is  entirely  elastic  so  £r  =  e(re)  and  ^  =  0  f** 
h(r)  =  0  Using  the  same  argument  for  the  other  strain  components  leads  to 

(4.13a) 
(4.13b) 
(4.13c) 


er=e(re)  +  e‘p) 

£e  =  £ee>  +  eeP) 


e  =ele)  +  e'zp) 


The  plane  strain  condition  dictates  that  ez  =  0,  so  using  Equations 
the  plastic  and  total  strains  are  given  by 
e<p)  =  0 
M^p>  =  e(ze) 


4.11  and  4.13, 


(4.14a) 

(4.14b) 

(4.14c) 


and 


e,  =  e?» 


c  —  4 _ 

ee  “  ee  +  z 


ez  =  0 


(4.15a) 

(4.15b) 

(4.15c) 


Hooke's  law  relating  the  elastic  strains  to  the  stresses  is  given  by 

Ee(re)  =ar-v(a8  +  az) 

E  4e)=ce-v(cz  +  ar) 

Eeze)  =  crz  -v(ar  +  ae) 


(4.16a) 

(4.16b) 

(4.16c) 


Eliminating  o.  from  Equations  4.16  via  the  yield  condition,  and  substitution  of  the 
resulting  elastic  strains  into  Equations  4.15  gives  the  following  expressions  for  the 
radial  and  circumferential  stress  components  in  terms  of  the  radial  and  circum- 
ferential  strain  components. 

[MN  + 1-  (M  +  1)(N  +  l)v]o,  =  [MN  + 1-  (M + N)v]2Ger 

+M(N  +  l)v2Gee  +  (M-l)vc„  W.17; 
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(4.17b) 


[MN  + 1-  (M  +  1)(N  + 1)  v]ae  =N(M  +  l)v2Ger 

+MN2Ge0  -[l-(M  +  l)v]au 


As  discussed  in  Section  3,  the  equation  of  equilibrium  is  given  by 


r^  +  Gr-ce  =  0 

dr 

(4.18) 

and  the  strain-displacement  equations  are 

du 

(4.19a) 

£r"  dr 

u 

e"  =  7 

(4.19b) 

Substitution  of  the  stress  Equations  4.17  into  the  equilibrium  Equation  4.18  and 
using  the  strain-displacement  Equations  4.19,  leads  to  the  following  differential 

equation  in  terms  of  the  radial  displacement  u, 

„d2u  du|\  (M-N)v  ,1  „2  P2  (1-2v)o 

rd?  +  r^L1+  MN  M'PU  MN  2G 

(4.20) 

where 

MN 

P  MN+1-(M  +  N)v 

(4.21) 

The  solution  to  Equation  4.20  is 

2Gu  =  A^1  +  A2r~Yi!  +  Ar 

(4.22) 

where 

A”-1-2Nv°u  =-(1-2v)Pb 

(4.23) 

(N-M)v  2  l(N-M)V  +  1 

7l  ^  2MN  P  ]j  (2MN)2  P2 

(4.24) 

(M-N)v  02  (M-N)V+J_ 

72  P  2MN  P  \  (2MN)2  p2 

(4.25) 

The  strains  can  be  computed  from  the  displacement  Equation  4.22  and  the  strain- 
displacement  Equations  4.19 
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2Ger  =  yA^1'1  -  Y2A2r'Yi_1  +  A 
2Gee  =  A^1'1  +  A2r'Y2_1  +  A 


(4.26a) 

(4.26b) 


Substituting  the  strain  Equations  4.26  into  the  stress  Equations  4.17  yields 
cr  =  C^A^1'1  +  cr2A2r~Y2“1  -Pb 
ae  =  CeiA1rY‘-1+Ce2A2r^-1-pb 


(4.27a) 

(4.27b) 


where 


Crl  =  -{[MN  + 1-  (M  +  N)v]Yi  +  M(N  +  l)v} 

c 

Cr2  =  I  {-[MN  + 1  -  (M + N)  v] y2  +  M(N  + 1) v} 
Cei  =  -i[N(M  +  1Wi  +  MN] 

C  =  MN  +  1-(M+  1)(N  +  l)v 


(4.28a) 

(4.28b) 

(4.28c) 

(4.28d) 

(4.28e) 


At  this  point,  the  unknowns  are  the  two  constants,  Ax  and  A2,  and  the  elastic- 
plastic  interface  radius  R.  Two  of  the  conditions  used  to  evaluate  the  unknowns  are 
provided  by  enforcing  continuity  of  radial  stress  and  radial  displacement  (or, 
equivalently,  circumferential  strain)  at  R.  The  elastic  zone  stresses  are  given  in  the 

elastic  solution  presented  in  Section  2.3,  with  the  substitutions  a  =  R  and 
pa  =  -cr(R).  Applying  the  yield  condition  of  Equation  4.2  at  r  equals  R  gives 

(A  oq 

or  (R)  =  -2(1-  Nv)pb  +  au 

The  circumferential  strain  at  r  equals  R  is  given  by 

2Gee  (R)  =  -2(N  - 1)  vpb  -  au  (4‘3C 

Using  Equations  4.26,  4.27,  4.29  and  4.30,  the  continuity  conditions  are  then  given 


+  Cr2A2R"Y2_1  =  2(l-Nv)Apb 
AxRYl_1  +  A2RY2_1  =  2(N  -  l)vApb 


(4.31a) 

(4.31b) 


where 

Apb  =Pb-Pb 


(4.32) 
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Equations  4.31  constitute  a  system  of  two  equations  in  terms  of  two  unknowns 
A^1-1  and  A2R‘T2_1. 

f  AjR71"1 1  =  j2(l-Nv)Apbl  (4  33) 

I  AoR'Y2_1 1  [2(N-  l)vAPb  J 


CrX  cr2 


1  1 

which  can  be  inverted  to  give 
fBil  _  f  AtRYl'1 1 


Bo 


A9R'Y2_1I  Crl-Cr2l_-1  crl  J 


1  -c 


r2 


(2(l-Nv)Apbl 
1 2(N-  l)vApb  J 


(4.34) 


where  Bx  and  B2  are  constants. 


The  elastic-plastic  interface  radius  R  is  found  by  enforcing  the  pressure  boundary 
condition  at  the  edge  of  the  hole 

,  x  (4.35) 

ar(a)  =  -pa 

which  leads  to  the  nonlinear  equation 


(  a  Y,_1 

CrlBl(  | J 


+  Cr2B2|  r 


-Y2-1 


~Pb  =“Pa 


(4.36) 


which  can  be  solved  for  the  radius  ratio  (a/R),  which  in  turn  gives  R.  Once  R  is 
determined,  the  constants  Ax  and  A2  can  be  recovered  from  B,  and  B2,  respectively. 


The  solution  presented  above  is  valid  with  continued  loading  as  long  as  or  <  az  in 
the  plastic  zone.  Florence  and  Schwer  performed  numerical  studies  of  the  stress 
difference  az  -ar  and  found  that  ctz  -ar  decreases  as  pb  increases,  and  becomes 
zero  first  at  the  edge  of  the  hole.  We  now  derive  the  applied  pressure,  pb,  and 
elastic-plastic  radius,  R,  associated  with  zero  stress  difference  at  the  edge  of  the 
hole. 


Substituting  Equation  4.34  into  Equations  4.27  gives 
,  r  vTi-1 

a.  =  CrlB, 


+  Cr2B2 


f  r  V*"1 


U ) 


R 


n-i 


=  CelB1^  +Ce2B2l  —  J  -Pb 


-Pb 

-Y2-1 


The  out-of-plane  stress,  az,  can  be  determined  from  the  yield  condition  and 
Equations  4.37. 


(4.37a) 

(4.37b) 
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_  _  ^ei  g  ( _L 

z"  N  BlU 


+  nB2IrJ  ~2vPb 


We  now  introduce  two  new  constants,  Bx  and  B2,  defined  by 
Bx  =  BxApb 
B2  =  B2Apb 

or  from  Equation  4.34 

ftl.  f1  -c'!ll2(1;NyH 


|bJ  Crl-Cr2  [-I  CrlJ\2(N-l)vJ 

Substituting  Equations  4.39  into  Equations  4.37a  and  4.38  gives 


ar=CrlBxApb-  +  Cr2B2Apb I  —  I  -Pb 


(4.38) 


(4.39a) 

(4.39b) 


(4.40) 


(4.41a) 


-2vp 


(4.41b) 


We  denote  as  p'b  and  R',  the  applied  pressure  and  radius  to  the  elastic-plastic 
interface,  respectively,  when  or(a)  =  oB(a).  Using  Equations  4.41,  and  setting 
ar(a)  =  az(a)  =  -pa  gives 


CrlBlAp'b^J  +CrM>{K)  =Pb_Pa 

.(  a  V1'1  .  cfl2  A  A_,(  a  V7*"1  - 9vA 


=2vp‘'p 


where 


APb  =  Pb  “  Pb 

Equations  4.42  constitute  a  system  of  two  equations  in  terms  of  unknowns 
Ap^y1  1  and  Ap'^)  ^  which  can  be  inverted  to  give 


__lf  ce2B2  -NCr2B2  |  Pb-Pa 
-1_  D  -CexBx  NCrlBx  Jl2vpb-pt 


(4.42a) 


(4.42b) 


(4.42c) 


^It7 

|Ap'bf#7 


(4.43a) 


where 


D  =  BxB2(CrXCe2-Cr2Cex) 

Equation  4.43a  can  be  rearranged  to  give 


(4.43b) 
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(4.44a) 


a  Y*'1 

kR^J  _  B1[Ce,(pb-p.)-NCI,(2vpb-P.) 
'  o  V"1  B2[Ce2(pb-p.)-NC,2(2vpb-p.) 


ap; 

'  a  N 

,R 7J 

1 

-d 

M 

1 

ap; 

a  " 

U', 

rH 
!  1 

3^ _ 

or 


<R/  \Yl+Y2 
,  a  j 


Bx 


B, 


Cei(Pb-Pa)-NCrl(2VPb-Pa). 

Ce2(Pb-Pa)-NCr2(2Vpb-pa) 


Once  —  is  determined,  Apb  can  be  obtained  from  Equation  4.43a. 
a 

,  _  (C92  ~  2  vNCr2  )pb  —  (Ce2  —  NCr2)pa  R  ^ 

Pb  "  B1(CrlCe2-CraC.l)  V  a  J 
Using  Equation  4.42c,  p'b  can  be  determined  from  Ap;. 


(4.44b) 


(4.45) 


When  the  free  field  yields,  R'  goes  to  infinity  implying  that  the  denominator  of  the 
right  side  of  Equation  4.44b  must  approach  zero.  Note  that  for  Nv  <  -,  the  free  field 

will  not  yield.  Numerical  studies  indicated  that  the  term  B2  is  nonzero  for  Case  II 
material  parameters.  This  implies  that  the  term  in  the  denominator  of  the  right 

side  of  Equation  4.44b  which  can  go  to  zero  (for  Nv  <  — )  is 

c.2  (pb  -  P. )  -  NCi2  (2vp„  -  p, )  =  0  (4-46) 

Further  numerical  studies  indicated  that  when 


f-  2NvCr2  Ce2  _  «■* 
P*  Pb  NCr2-Ce2  Pb 


(4.47) 


a  solution  for  R'  using  Equation  4.44b  is  possible;  thus  yielding  of  the  free  field 

occurs  after  the  formation  of  the  second  plastic  zone.  Therefore,  the  condition 

pa  <  p*b  defines  Case  Ha.  Conversely,  when  pa  >p*b,a  solution  for  R'  using  Equation 

4.44b  is  not  possible,  indicating  that  yielding  occurs  prior  to  the  formation  of  the 
second  plastic  zone.  Therefore,  the  condition  pa  >  pb  defines  Case  lib. 


For  Case  lib,  the  solution  presented  above  can  be  modified  to  include  the  case  of  pb 
equal  to  pb.  Arguments  similar  to  those  presented  in  Section  3.4  may  be  used  to 
show  that  as  pb  approaches  pb,  the  constant  Ax  in  Equations  4.22,  4.26,  and  4.27 

must  go  to  zero  and  that  R  must  approach  infinity. 
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Thus  when  pb  equals  pb,  the  radial  stress  is  given  by 

=Cr2A2r_Y2'1-Pb 

The  constant  A2  can  be  determined  from  the  pressure  boundary 
edge  of  the  hole,  crr(a)  =  ~pa 
Cr2A2a”1f2”1  _Pb  =  “Pa 


(4.48) 

condition  at  the 

(4.49a) 


or 


A  - 

a2-c 


-Y2-1 


yr2e 


(4.49b) 


The  remainder  of  the  solution  presented  above  for  the  plastic  zone  is  valid  for  pb 
equals  pb. 


4.2  CASE  Ha,  PHASE  3  SOLUTION. 


In  this  section,  we  present  the  solution  for  Phase  3  of  Case  Ila  which  contains  two 
plastic  zones  and  an  outer  elastic  zone  (Figure  4-3).  We  also  examine  the  case  where 
pb  is  exactly  equal  to  pb  and  the  elastic-plastic  boundary  goes  to  infinity. 


Florence  and  Schwer  postulated  that  with  increased  loading  (p„  >p’h)  a  radius 
R  >  a  exists  where 

ar(R)-aI(R)  =  0 


In  the  outer  plastic  zone,  we  still  have  a,<a,<a„  the  same  condition  as  in  the 
Phase  2  solution  presented  in  Section  4.1.  In  fact,  the  stress  and  strain  fields 
developed  in  Section  4.1  still  apply.  For  convenience,  they  are  repeated  be  ow 


(Equations  4.26  and  4.27). 

2G€r=Y1A1rv‘-1-72A2r-T2-1  +  A 

2Gee  =  AxrTl_1  +  A2r'T2_1  +  A 
ar  =  C^A^1"1  +  Cr2A2r"V2_1  -  pb 
ce  =  CelAjrYl-1  +  Ce2A2r-Y2-x  -  pb 


(4.26a) 

(4.26b) 

(4.27a) 

(4.27b) 
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The  continuity  conditions  at  the  elastic-plastic  interface  R  developed  in  Section  4.1 
still  apply.  These  conditions  led  to 


[Bi\  _  J  AxRYl_1 1 


Bo 


A2R_Y2_1 


1  -c 


r2 


Crl-Cr2  l-l  CrlJ 


[2(l-Nv)Apbl 
i  2(N  -  l)vApb  J 


(4.34) 


However,  the  elastic-plastic  interface  R  is  no  longer  found  by  enforcing  the  pressure 
boundary  condition  at  the  edge  of  the  hole  (Equations  4.35  and  4.36).  Instead,  we 
use  the  condition  of  Equation  4.50.  The  out-of-plane  stress  cz  is  determined  (in 
terms  of  oe  and  Gu)  from  the  yield  condition  of  Equation  4.2,  and  Equation  4.27b  is 
substituted  for  oe  to  yield  the  following  nonlinear  equation 


r  b 
.  rl  N  H 


^rV’-1 


+  C--W 


VRy 


-(l-2v)pb  =  0 


(4.51) 


which  can  be  solved  for  the  radius  ratio  (R/R) .  The  solution  in  the  inner  plastic 
is  required  to  define  the  individual  values  of  R  and  R. 


zone 


In  the  inner  plastic  zone,  Florence  and  Schwer  postulated  that  the  stress  condition 
is  oe  <  ar  =  az  so  that  yielding  is  governed  by 


fx  =  ae-Nor+ou  =0 
f2  =  ae-Noz+ou  =0 


(4.3a) 

(4.3b) 


From  Equation  4.3a, 

o6  =  Nor-ou 

so  the  equilibrium  equation  (Equation  4.18)  can  be  written  as 


(4.52) 


r^-(N-l)Gr  +  ou=0 
dr 

The  solution  to  Equation  4.53,  with  Gr(a)  =  —  pa,  is  given  by 


(4.53) 


N- 1 


(  r  V*-1 


Va  ; 


N-l 


(4.54) 


Continuity  of  the  radial  stress  at  the  plastic-plastic  interface  (r  equal  to  R)  leads  to 


CA 


'rV'-1 

vRy 


+  Cr2B2 


^vTrl 

vRy 


-Pb=-[Pa  +  NnJ 


Rf-1 


N-l 


(4.55a) 


or 
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EY'"1 


sr-1 


N(l-2v)pb-(N-l)C„B1|-|  -(N-1)Ci2B2^rJ 

-J  -  (N-l)pa+ou 

which  can  be  solved  for  R  since  Bt>  B2,  and  (R/R)  have  been  determined  previously. 
Once  R  is  determined,  R  is  computed  from  (R/R)  and  the  constants  A,  and  A,  can 
be  recovered  from  B,  and  B„  respectively.  Thus  the  solution  in  the  outer  plastic 
zone  is  complete. 


In  the  inner  plastic  zone,  the  strain  rates  are  given  by 

e,pl  =  + X2  =  -M>.,  (4'56a) 

r  dct  oar 

iv  =  x*ii  +  K‘Si  =  -MX2  (4'56C) 

dcrz  daz 

From  the  flow  rule  of  Equation  4.56,  it  follows  that 

e<p)+Me<p)  +  e(zp)=0  (4,57) 

which,  when  integrated  with  respect  to  time,  gives 

E(p)  +  M£(ep)  +  e(2p)  =h(r)  (4.5S) 

Let  the  plastic  strains  be  e(rp) ,  e<p) ,  and  elP)  at  radius  r  when  the  load  pb  has  placed 

the  plastic-plastic  radius  at  r.  Then,  when  R  equals  r, 

e<p)  +  Me<p)  +  ezp)  =  h(r)  (4,59) 

In  the  outer  plastic  zone,  the  plastic  strain  relations  of  Equation  4.11  still  hold  so 
h(r)  must  be  zero.  Therefore  Equation  4.58  becomes 

e(p)  +  Me‘p)  +  e(zp)  =0  (4’60) 

Snhstitutinff  Eauation  4.60  into  the  plane  strain  condition  (ez  =  e(ze)  +e(2p)  =  0)  gives 


(4.56b) 


(4.56c) 


(4.57) 


(4.58) 


(4.59) 


Substituting  Equation  4.60  into  tne  plane  swam  tuauiuu* 
er  =  -Mee  +  Me(ee)  +  e(re)  +  e(2e) 

Using  Equation  4.52,  and  or  =  o2,  reduces  Hooke's  law  to 
Ee(re)  =  Ee(ze)  =  [l-  (N  +  l)v]ar  +  vou 
Ee^e)  =(N-2v)ar-au 


(4.61) 


(4.62a) 

(4.62b) 
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where  cr  is  given  by  Equation  4.54.  Substitution  of  the  elastic  strains  of  Equation 
4.62  into  Equation  4.61,  followed  by  the  substitution  of  the  resulting  er  into  the 

compatibility  equation 


r_i  +  e0-er  =0 
dr 


leads  to  the  equation 


2G£( 


rM+le0 


'  MN  +  2-2(M  +  N  +  l)vr,„ 

(N-l)(l+v)  l' 


M+N-l 


if 

°uj  QN-1 


(4.63) 


(M  +  2)(l-2v)au 


(4.64) 


(N-l)(l+v) 

The  differential  Equation  4.64  is  directly  integrate.  Integration  from  r  to  R  gives 


2Gee  =  2Gee(R) 


-\M+1 


MN  +  2-2(M  +  N  +  1)v[(N_1)  ] 

(N-l)(l  +  v)  11  1 


-xM+N 

*  -i 


N-l 


(M  +  2)(l-2v)ou 


r- 


(4.65) 


(N-  1)(M+ 1)(1+ v)| 

where  2Ge0(R)  can  be  determined  from  Equation  4.26a.  Substitution  of  the  elastic 
strains  from  Equation  4.62  into  Equation  4.61  results  in 


2Gf  =  -2GMe0  + 


MN  +  2-2(M  +  N  +  l)v 


1+  v 


M-2v 
1+ v  G" 


(4.66) 


where  or  is  given  by  Equation  4.54  and  2Ge0  is  given  by  Equation  4.65. 


When  the  increasing  pb  approaches  pb  (the  load  which  initiates  free-field  yielding), 
we  can  use  the  same  arguments  as  presented  in  Section  3.4  to  show  that  the 
constant  A1  in  Equations  4.22,  4.26,  and  4.27  must  go  to  zero,  and  that  R  must 

approach  infinity. 

Thus,  when  pb  approaches  pb,  our  unknowns  are  reduced  to  A2  and  R.  The  con¬ 
dition  of  Equation  4.50,  with  oz  eliminated  via  the  yield  condition,  and  Equation 
4.27b  substituted  for  oe  (with  Aj  equal  to  zero),  gives 
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(cr2  -  )a2R-’'1-'  -  (1-  2v)pb  =  0 

which  can  be  solved  for  A.R"'--’.  Continuity  of  radial  stress  at  R  leads  to 

_  ^  r  a  IfRf'1  ,  qu 

Cr2A2R"Y2  1  “ Pb  =  “I  Pa  +  N _ i  \  a  J  N-l 


(4.67) 


(4.68) 


N-l 


which  can  be  solved  for  R  since  A2R—  was  determined  Previously^  Once  the 
unknowns  are  determined,  the  stress  and  strain  expressions  developed  p„  P„ 

are  still  valid  for  Pb  =  Pb* 

4.3  CASE  lib,  PHASE  3  SOLUTION. 

The  solution  for  Phase  3  of  Case  lib  is  presented  in  this  section.  As  shown  in  Figure 
4-4,  in  this  phase  we  have  an  inner  plastic  zone  with  conditions  e 

governed  by  the  yield  condition  ^  ^ 

and  a  to  roter^fc'zone  with  conditions  <J,  =  o.  <  governed  by  the  yield  con- 


ditions 


f  =ar-Naz  +  au  =  0 

A  (4 

f2  =  ae-Naz+au=0 

The  conditions  in  the  far  outer  plastic  zone  were  postulated  to  exist  based  on 
numerical  studies.  From  the  external  pressure  boundaiy  condition  and  the  yie 
condition  of  Equation  4.5,  we  have 


(4.5a) 


(4.5b) 


(4.69) 


ar  =  c6  - -pb 

In  the  inner  plastic  zone,  the  stress  and  strain  fields  developed  in  Section  4.1  still 
apply.  The  equations  representing  these  fields  (Equations  4.26  and  4.27)  are  again 

repeated  for  convenience. 

-1  a  -7,-1  a  (4.26j 

2Ger  =  YiAir7l_1  -  Y2A2r  +A 

.  (4.261 

2Gee  =  A1rYl'1  +  A2r"Y2  +A 

(4  27i 

or  =  C^Ap"71-1  +  Cr2A2r'72  1  -  Pb 

-  .  _7„_x  .  (4.271 

ce  =  C61A1rYl_1  +  Ce2A2r  72  -pb 


(4.26a) 

(4.26b) 

(4.27a) 

(4.27b) 
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where  A  =  -(l-2v)pb. 


At  this  point,  the  unknowns  in  the  inner  plastic  zone  are  the  constants  Ax  and  A2, 
and  the  radius  of  the  boundary  between  the  inner  and  far  outer  plastic  zones, 
denoted  as  R.  Based  on  the  yield  conditions  of  Equations  4.2  and  4.5,  along  with 
Equation  4.67,  it  is  reasonable  to  postulate  that  c6(r)  =  -pb  in  the  inner  plastic 

zone.  This  condition,  along  with  continuity  of  ar  at  R,  results  in  two  equations  m 
the  two  unknowns  A.R*1'1  and  A2R'Y2'\  which  can  be  solved  to  yield 


Yi-l 


A2R"Y2_1 


crlce2-cr2c 


r2^ei 


J62 


— c 


ei 


where  B1  and  B2  are  constants.  After  solving  for 
of  Equation  4.27  can  be  rewritten  as 


\Yi-l 


or  =  CrlBx 


R. 


+  Cr2B2 


r 

,R 


\-Y*-l 


"Pb 


Ge  —  (-'ei® 


xIrJ 


+  Ce2B2|  y 


-Y2-l 


"Pb 


-c 


r2 


'rl  J 


Apb 

lAPb. 


(4.70) 


Bx  and  B2,  the  stress  expressions 


(4.71a) 

(4.71b) 


At  the  edge  of  the  hole,  we  have  the  pressure  boundary  condition  oT(  a)  -  -pa  which, 
when  substituted  into  Equation  4.71a,  leads  to 

\-Y2-l 


Gr(a)  =  CrlBx 


f  a  V‘_1 

i  +CA 


a_ 

R. 


-Pb  =“Pa 


Equation  4.72  is  a  nonhnear  equation  which  can  be  solved  for 


a 

R. 


(4.72) 

which  in  turn 


yields  R.  Once  R  is  determined,  the  constants  Ax  and  A2  can  be  recovered  from  Bx 
and  B2.  The  stress  and  strain  fields  in  the  inner  plastic  zone  are  now  completely 

determined. 


In  the  far  outer  plastic  zone,  the  stress  and  strain  fields  developed  in  Section 
apply.  The  stress  fields  are  given  by 

Gr  =  “Pb 
Ge  =  "Pb 

Gz=^K-pb) 


3.4 

(4.73a) 

(4.73b) 

(4.73c) 
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It  was  shown  in  Section  3.4  that  the  sum  of  the  radial  and  circumferential  plastic 

A 

strains  was  a  constant,  C ,  given  by 


£(p)  +  £flP) 


1  r  Pb+£iL  +  2vp1=°--<1-2N%  =  C  (4.74) 

‘EmL'N  N  PlJ  2GMN(l+v) 

The  equation  of  compatibihty,  strain  decomposition,  and  Hooke's  law  lead  to  the 
following  differential  equation  for  £(ep) 


rMUe!.p)-£(p)  =0 

dr 


which,  using  Equation  4.74,  can  be  solved  to  give 


elp)  =  — + A, 


.-2 


(4.75) 


(4.76) 


where  A3  is  a  constant.  It  was  also  shown  in  Section  3.4  that  the  postulation 
ae(R)  =  -pb,  Hooke's  law,  and  continuity  of  total  circumferential  strain  require  that 


„<p) 


must  be  continuous  across  R .  This  condition  is  used  to  determine  A3 . 

4p>(r)=|+a3r 


or 


A,  =  R5 


£(p)_C 
£e  2 


(4.77a) 


(4.77b) 


where  e^(r)  is  computed  from  the  solution  in  the  inner  plastic  zone.  Equation 

4.26b,  evaluated  at  R,  gives  the  total  strain  ee(R).  Equations  4.27  and  4.2  can  be 

used  to  compute  the  stresses  at  R.  Hooke’s  law  can  then  be  applied  to  compute  the 
elastic  strain  e^R).  Finally,  the  plastic  strain  £<p)(r)  is  computed  by  subtracting 

£(ee)(R)  from  £0(r). 


The  solution  developed  above  is  valid  until  cz  equals  Gr  at  r  equals  a.  When  this 
occurs,  a  third  plastic  zone  is  formed  which  begins  at  the  edge  of  the  hole  and 
migrates  outward.  When  the  third  plastic  zone  is  present,  the  Phase  4  solution 
(which  will  be  presented  in  Section  4.4)  applies.  However,  since  the  Phase  4  solution 
applies  to  both  Cases  Ha  and  lib,  we  present  here  the  derivation  for  the  load  pb 
which  will  cause  the  formation  of  a  third  plastic  zone  for  Case  lib.  The  load  pb  and 
plastic-plastic  boundary  radius  R  when  the  third  plastic  zone  begins  to  form  will  be 
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denoted  as  p"  and  R" ,  respectively,  and  are  derived  below.  The  derivation  is 
similar  to  that  given  in  Section  4.1  for  p'b  and  R'. 


As  before,  we  introduce  two  constants  Bj  and  B2,  defined  by 
Bx  =  BxAPb 
B2  =  B2Apb 

From  Equations  4.70  and  4.78, 


[b,1 

1  1  J 

f  c,„  -cr2 1 

Uj 

[  CrlCe2  —  Cr2Cei  1 

[-C.1+C.J 

Equation  4.71 


r>  (  t*  Yfl_1  C 
„  i  +  — B.i  — 

•  N  HrJ  n  "U. 


-Y2-1 


-2vpb 


Substituting  Equation  4.79  into  Equations  4.71a  and  4.80  gives 


( r  V1"1  -  ( r 

ar  =  C^BbApbjjrJ  +  Cr2B2APb|j- 


-Y2-1 


-Pb 


f  r  V*2'1 


u 


-2vpb 


Using  Equations  4.81  and  setting  ar(a)  =  cz(a)  =  -pa  gives 

\Yi_1 


C^B^p" 


a 


C 

N 


R" 

£ 

R" 


+  Cr2B2Ap 


f  a  VY2_1 


vFj  =Pb'p- 


a  1Y'-1  C 


^b.ak  £  b2ap;1  tz 


_a 

,R" 


\~12~1 


— 2vpb  Pa 


where 

APb=pb-Pb 

Equations  4.82  constitute  a  system  of  two  equations  in  terms  of  unknowns 
and  APb^J  which  can  inverted  to  give 


\Yi-l 


Ap<f) 


_Y2-1 


Ce2B2  -NCr2B 


22^2 


D  [-CejBj  NCrlBx  Jl2vPb-P 


Pb-Pa 


(4.78a) 

(4.78b) 

(4.79) 

(4.80) 

(4.81a) 

(4.81b) 

(4.82a) 

(4.82b) 

(4.82c) 

(4.83a) 
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where 


(4.83b) 


D  =  B1B2(CrlC62-Cr2Cei) 

Equation  4.83a  can  be  rearranged  to  give 


f  o  V**-1 

> 

\^Z± 

u 

Jr, 

(  a  V1'1 

Ap'b 

(jr) 

®i 


B, 


C01(ph-pa)-NCrl(2vpb-pj; 
Ce2  (pb  -  Pa  )  -  NCr2  (2  vpb  -  Pa  ) 


(4.84a) 


or 


R'-Y'*1''  B,[C„(pt-P.)-NCrl(2vpt-p.)] 


Once 


R 

V  a  J 
fR'A 


\  a  j 


IS 


B2[CM(pb-p.)-NCf2(2vpb-p.)] 
determined,  Apb  can  be  obtained  from  Equation  4.83a 


(4.84b) 


(Cfl9  -  2  vNCr2  )pb  -  (C92-NCr2)Pa 
APb_  BjCA-CrAi) 


rR"Vl_1 


v  a  J 


(4.85) 


and  p^  can  be  determined  from  Equation  4.82c. 

4.4  CASE  Ha  AND  lib,  PHASE  4  SOLUTION. 

In  this  section  we  present  the  solution  for  Phase  4,  which  contains  three  plastic 

(Figure  4-5).  Again,  this  solution  is  only  applicable  for  Nv  <  ^  d*scusse^ 


zones 

earlier,  yielding  begins  in  Case  Ua  at  the  edge  of  the  hole  with  an  inner  plastic  zone. 
With  increased  loading,  a  second  plastic  zone  forms  which  also  migrates  outward 
from  the  edge  of  the  hole.  The  initial  plastic  zone  is  now  the  outer  plastic  zone,  whi 
the  outermost  region  is  elastic.  In  Phase  4,  the  outermost  region  becomes  plastic 
when  pb  becomes  large  enough  to  cause  the  free  field  to  yield. 

In  Case  lib,  yielding  begins  at  the  edge  of  the  hole  with  an  inner  plastic  ™ne,  with 
the  outermost  region  being  elastic.  This  is  followed  by  free-field  yielding.  In  Phase 
4  a  third  plastic  zone  is  formed  which  begins  at  the  edge  of  the  hole  and  migrates 
outward.  So,  although  the  path  taken  to  get  to  Phase  4  differs  in  Cases  Ha  and  I  , 
the  number  of  plastic  zones  and  the  conditions  in  each  zone  are  identical,  and 

the  solution  is  identical  for  both. 
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The  stress  condition  in  the  inner  plastic  zone 
governed  by 

f\  =  ae-Nar+au  =0 
f2=Ge-NGz+Ga=0 

The  stress  condition  in  the  outer  plastic  zone 
condition  is 


(a  <  r  <  R)is  Ge  <  ar  =  gz  with  yielding 

(4.3a) 

(4.3b) 

(1  <  r  <  R)  is  Ge  <  Gr  <  gz  so  the  yield 


L  _ 

Finally,  in  the  far  outer  plastic  zone  (R  <  r  <  °°),  we  postulate  that  the  stress  con¬ 
dition  is  Ge  <  Gr  =  Gz  so  yielding  is  given  by 


fx  =  Gr-NGz  +  Gu  =0 
f2  =  g6-Ngz+gu  =0 

The  external  pressure  boundary  condition  and  the  yield 
lead  to 


(4.5a) 

(4.5b) 

condition  of  Equation  4.5 


<*r  =  <*0  =  "Pb 


(4.69) 


As  discussed  in  sections  4.1  and  4.2,  the  stress  and  strain  fields  in  the  outer  plastic 
zone  are  give  by 

2Ger  =  YiA^1"1  -  Y2A2r"Y2"1  +  A 

(4.26a) 

2Gee  =  Ap11'1  +  A2r~Y2_1  +  A 

(4.26b) 

Gr  =  CrlA1rYl-1  +  Cr2A2r"Y2_1  -  pb 

(4.27a) 

g0  =  CelA1rYl_1  +  Ce2A2r_Y2_1  -  pb 
where  A  =  -(l-2v)pb. 

(4.27b) 

At  this  point  the  unknowns  are  the  constants  Ax  and  A2  and  the  plastic  zone 
boundary  radii  are  R  and  R.  Just  as  in  Sections  3.4  and  4.3,  we  postulate  that 
g6(r)  =  _Pb  in  the  outer  plastic  zone.  This  condition,  along  with  continuity  of  Gr  at 

R ,  results  in  two  equations  for  two  unknowns  A^1'1  and  A*RT»“\  which  can  be 
solved  to  yield 


b2 


j  1 

1a2b-"-xJ 


^rl^e2  —  ^r2^01 


(4.70) 
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where  Bx  and  B2  are  constants.  After  solving  for  Bx  and  B2  ,  the  stress  expressions 
of  Equation  4.27  can  be  rewritten  as 

'-72-1_£b  (4.71a) 


cJr  =  CrlBx 

=  ^ei®i 


(  r  V1'1 


+C'A  I. 


\-TT2"1 


-Pb 


(4.71b) 


In  the  inner  plastic  zone,  the  stress  and  strain  fields  developed  in  Section  4.2  apply. 
The  radial  stress  in  that  zone  is  given  by 

(4.54) 


*r-  ya  +  ^_x 


nN-1  /t 

i) 

aJ  N-l 


(4.50) 


In  Section  4.2,  we  also  postulated  that  at  the  inner  boundary  of  the  outer  plastic 
zone  (R) 

ar(R)-az(R)  =  0 

Using  Equation  4.50  and  the  yield  condition  Equation  4.2  gives 

Ncr  (R)  -  o0  (R)  -  ou  =  0  (4,86) 

Substituting  Equation  4.71  into  Equation  4.86  gives  the  following  nonlinear 
expression 

1  ✓ -  \  -Vo-1 

(4.87) 


^RVrl 

R 


+  (NC,2-C92)B2[|]  -N(l— 2v)pb  =0 


(NC^-C.OB^ 

which  can  be  solved  for  the  radius  ratio  (r/r).  Using  Equations  4.54  and  4.71a  to 


enforce  continuity  of  or  at  R  gives 

f  tT 


rTY2_1 


or 


_,N-r  N(l-  2v)pb  —  (N— l)CrlB, 


R 

a 


^Rv‘ 
\R  j 


O.  1 

+  °u 

N-lJ\ 

,aj 

N-l 

L-l 

-  \  y* 

v,  rin 

(4.88a) 


■  (N  -  l)Cr2B2 


R 


(4.88b) 


(N-l)pa+cu 
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which  can  be  solved  for  R  since  B„  B2>  and  |  have  been  determined  previously. 

After  solving  for  R,  R  can  be  recovered  from  £  and  A,  and  A2  can  be  recovered 
from  Bx  and  B2  respectively. 

The  solution  in  the  far  outer  plastic  zone  is  identical  to  that  given  in  Section  4.3, 
Equations  4.74  through  4.77. 
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SECTION  5 

PROCEDURES  FOR  APPLYING  EXTENDED  SOLUTIONS 

In  this  section,  the  procedures  for  applying  the  extended  solutions  are  summarized 
in  13  tables.  These  tables  are  described  briefly  in  the  following  paragrap  s. 

Because  of  the  large  number  of  multi-page  tables  and  the  relatively  sma  amoun  o 
accompanying  text,  all  tables  are  placed  at  the  end  of  the  section. 

The  first  step  in  the  solution  process  is  to  define  the  problem.  The  material  an<* 
load  parameters  that  must  be  known  prior  to  solving  the  problem  are  listed  m  Tab 
5-1.  Restrictions  on  the  known  parameters  (discussed  in  Section  2)  are  a  so 

included  in  that  table. 

The  second  step  is  to  determine  which  subcase  (Case  la,  lb,  Ila  or  lib)  applies.  The 
Case  I  and  Case  II  definitions  were  developed  in  Section  2.  The  definitions  o  ases 
la  and  lb  were  developed  in  Section  3.1,  and  the  definitions  of  Cases  Ila  and  lib 
were  developed  in  Section  4.1.  These  definitions  are  summarized  m  Table  5-  . 

Using  only  the  parameters  listed  in  Table  5-1,  one  can  use  Table  5-2  to  determine 
whether  a  particular  problem  falls  into  Case  la,  Case  lb,  or  Case  II  However,  an 
additional  parameter  (pb)  is  required  to  distinguish  Case  Ila  from  Case  lib. 

A  list  of  "computed  parameters”  (which  can  be  calculated  from  the  known 
parameters  listed  in  Table  5-1  and  are  commonly  used  in  the  various  solutions)  is 
provided  in  Table  5-3.  Equation  numbers  from  the  earlier  derivations  are  listed. 
Where  the  same  equation  was  used  in  Sections  3  and  4,  both  equation  numbers  are 
shown.  Some  of  the  parameters  contained  in  Table  5-3  pertain  to  all  plastic 
solutions;  others  pertain  only  to  Case  II  (e.g.,  pb>;  while  still  others  pertain  only  to 

either  Case  Ila  of  lib.  We  have  found  it  useful  to  compute  and  save  these 
parameters  at  this  point  in  the  solution  process. 

Once  the  appropriate  subcase  (Case  la,  lb,  Ha,  or  lib)  is  determined  the  next  step  is 
to  determine  which  phase  of  the  solution  is  applicable.  As  discussed  at  the 
beginning  of  Section  3,  Case  la  consists  of  two  phases,  while  Case  lb  consists  of  up 
to  four  phases.  Both  Cases  Ila  and  lib  consist  of  up  to  four  phases  (Section  4).  We 
use  the  terminology  "up  to  four  phases"  since  in  the  cases  where  l-2Nn  <  0,  the  ee 
field  will  not  yield  and  only  the  first  three  of  the  four  phases  are  applicable.  F 
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cases,  the  applicable  phase  is  simply  a  function  of  the  value  of  the  external  pressure 
Pb- 

The  process  for  determining  the  proper  case  and  subcase  is  summarized  in  the 
flowchart  shown  in  Figure  5-1.  Each  phase  of  each  subcase  constitutes  a  branch  on 
the  flow  chart.  The  branch  symbols  used  in  the  flow  chart  are  defined  m  Table  5-4. 
Also  contained  in  Table  5-4  are  references  to  figures  in  Sections  3  and  4  that  show 
the  stress  conditions  in  each  region  for  each  branch. 


Once  the  subcase  and  the  value  of  the  load  pb  are  known,  one  can  use  Table  5-5  to 
determine  which  branch  of  the  flowchart  (i.e.,  which  phase  of  the  solution)  is 
applicable.  Tables  5-4  and  5-5  also  contain  references  to  a  "solution  table"  for  each 
branch.  The  solution  tables  (Tables  5-6  through  5-13)  contain  step-by-step  solution 
procedures  for  each  of  the  branches.  Thus,  once  the  branch  is  determined,  the 
solution  procedure  may  be  found  in  the  applicable  solution  table. 

The  format  of  the  solution  tables  is  as  follows.  First,  unknown  coefficients  and 
plastic  boundary  radii  are  computed  where  required.  Next,  stress  and  strain 
response  quantities  are  given  for  each  region  in  the  solution.  Finally,  an  expression 
for  hole  closure  (referred  to  in  the  tables  as  "tunnel"  closure)  is  given.  If  the  only 
response  quantity  of  interest  is  closure,  one  may  skip  the  middle  part  of  each 
solution  table.  That  is,  one  can  compute  the  required  coefficients  and  boundary  radn 
and  then  proceed  directly  to  the  closure  calculation.  Equation  numbers  are  again 
provided  for  reference.  Primed  equation  numbers  are  used  to  indicate  that  the  form 
of  the  equation  has  been  modified  from  that  shown  in  the  derivation.  In  the  solution 
tables,  an  equation  denoted  with  multiple  numbers  is  an  algebraic  combination  of 
the  equations  bearing  those  numbers. 

Many  of  the  solutions  require  solving  a  nonlinear  equation  of  the  form 
C^x*1-1  +  C2x_1f2_1  +  C0  =  0 

where  C0,  Ci,  C2,  Yi  and  Y2  are  constants  and  x  is  the  unknown.  We  have 
successfully  used  both  linear  iteration  and  Newton-Raphson  techniques  to  solve 
nonlinear  equations  of  this  form.  The  linear  iteration  technique  is  summarized 

below. 
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Figure  5-1.  Flowchart  for  determining  solution  case. 


Figure  5-1.  Flowchart  for  determining  solution  case  (Continued). 
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Figure  5-1.  Flowchart  for  determining  solution  case  (Continued). 
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An  initial  guess  for  the  unknown  x  is  required.  For  most  problems  and  y2  are 
approximately  one,  so  a  reasonable  initial  guess  for  x  can  be  obtained  from 


C1((o)x)°  +  C2((o)x)"2+Co=0 


(5.2a) 


or 


(°), 


-a 


'c1+c0 


(5.2b) 


where  the  left  superscript  on  x  denotes  the  value  of  the  iteration  counter  k.  We  then 
successively  solve  for  ^k+1^x  using 


(k+1)x  = 


-q-c^x 


(k)vY,-l 


-Y2-I 


(5.3) 


until  acceptable  convergence  is  achieved.  Convergence  is  measured  by 

|(k+i)Y_(k) 


(k  +  l). 


<8 


(5.4) 


where  8  is  a  suitably  small  number  (e.g.,  1(H). 
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Table  5-1.  Known  parameters. 


Elastic  parameters: 

(n  ^ 

Poisson's  ratio  n  I  0  <  v  <  — 

Shear  modulus  G  (or  E  =  2G(l+v)) 

Strength  parameters 

Unconfined  compressive  strength  cu  (au  >  0) 

Friction  angle  <|>  (<|>  >  0) ,  or  N  =  (N  >  x) 

Dilatancy  angle  co  (0  <  co  <  <t>) ,  or  M  =  1_gi~  (1<M<N) 

Loading: 

Internal  pressure,  pa  ^Pa  >  0  for  Nv  >  - ,  or  0  <  pa  <  f°r  NV  < 

External  pressure,  pb  (pb  >  0) 

Geometry 

Interior  radius  a  (a=l  to  normalize  results  with  respect  to  interior 

radius). 
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iH  |  <N 


Table  5-2.  Case  la,  lb,  Ha,  and  lib  conditions. 


Propertv  Relations 

Internal  Pressure 

Case 

Nv  <  (N  +  l)v  <  1 

.  .  vau 

Pa  l-(N  +  l)v 

lb 

1  <  Nv  <  (N+l)  v 

Pa  >0 

lb 

Nv  <  1  <  (N+l)  v 

Pa  >0 

la 

—  <  Nv  <  (N  +  l)v  <  1 

2 

n  -  vau  i 

Pa  l-(N+l)v 

Ila 

0  <  Nv  <  -^ 

2 

vc“  <P  < 
l-(N  +  l)v  Pa  l-2Nv 

lib 

0  <  Nv  <  - 
2 

va„ 

-  7T7  "  TV-  <  Pa  <  Pb 

1-(N+  l)v 

Ila 

Note:  p*  is  a  function  of  several  variables.  See  Equation  4.47.  This 
equation  and  other  contributing  equations  are  summarized  in 
Table  5-3. 
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Table  5-3.  Computed  parameters. 


2_  MN 

M  MN+1-(M  +  N)V 

(3.39),  (4.21) 

(N-M)v  2  (N-M 

Yl  ^  2MN  ^  \{  2MN  J  p2 

(3.42),  (4.24) 

„2(M-N)v,„2 

Yz  ^  2MN  ^  \l  2MN  J  p2 

(3.43),  (4.25) 

a  °°  if  Nv  <  - 

(2.28b) 

Pb"l-2Nv  2 

Pb  =  ~  if  Nv  >  i 

A  =  -l-2Nv°“=-(1-2V)Pb 

(3.41),  (4.23) 

crl  =  -{[MN  + 1-  (M  +  N)v]y!  +  M(N  +  l)v} 

c 

(3.46a),  (4.28a) 

Cr2  =  £  {-[MN  + 1  -  (M  +  N )  v]  y2  +  M(N  + 1)  v} 

(3.46b),  (4.28b) 

Cei  =  ^[N(M  +  l)vYl  +  MN] 

(3.46c),  (4.28c) 

Ce2=^[-N(M+1)VY2+MN] 

(3.46d),  (4.28d) 

C  =  MN  + 1-  (M  +  1)(N  +  l)v 

(3.46e),  (4.28e) 

a  cj„-pb(l-2Nv) 

EMN 

Apb  =Pb-Pb 

(3.95),  (4.74) 

(3.50),  (4.32) 

Case  II  only: 

..  .  2NvCr,-C92  (4-47) 

Pb  Pb  NCr2-C62 
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Table  5-3.  Computed  parameters  (Continued). 


Case  Ha  only: 

iB,l- 

B. 


2  J  ~Tl 

^  a  J 


C,i-Cr2 
Bx 


1  -Cr2lj2(l-Nv)l 
L-l  CrlJ\2(N-l)vJ 

Cfl1(ph-pa)-NCrl(2vpb-pa) 


B2[Ce2(pb-pa)-NCr2(2vpb-pa) 
,  _  (Ce2  -  2 vNCr2  )pb  —  (C02  -  NCr2  )pa  f  R_ 

Pb "  B1(CrlC62l  ^  a 

Pb  =  Pb“APb 


\Yi-l 


Case  lib  only: 


B,  I 


C.2-C 


r2 


I  Bo  i  CrlCG2  -  Cr2Cel  [  C01  +  Crl , 


R"Yl+Y2_  B1[Cel(pb-pa)-NCrl(2vpb-pa) 


v  a  y 


Ap;'  = 

Pb  =  pb-AP;' 


B2  [c92  (pb  -  Pa )  -  NCr2  (2  vpb  -  Pa ) 

(Ce2  ~  2 vNCr2  )pb  —  (Ce2  —  NCr2  )pa  R  ^ 
B^CrtC,,)  ^a' 


Ti-1 


(4.40) 

(4.44b) 

(4.45) 

(4.42c) 


(4.79) 

(4.84b) 

(4.85) 

(4.82c) 
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Table  5-4.  Identification  of  flow  chart  solution  branches. 


Pqcd 

Phase 

Branch 

Symbol 

Solution 

Table 

Solution 

Conditions 

Figure 

Number 

VttOv 

I  (a  or  b) 

i 

1-1 

5-6 

Elastic  throughout 

3-1 

T  fa  or  b) 

2 

1-2 

5-7 

One  plastic  zone, 
free  field  elastic 

3-2 

Tb 

3 

Ib-3 

5-8 

Three  plastic  zones, 
free  field  elastic 

3-3 

1  u 

lb 

4 

Ib-4 

5-9 

Four  plastic  zones, 

(free  field  yielded) 

3-4 

II  fa  or  b) 

1 

II-l 

5-6 

Elastic  throughout 

4-1 

II  fa  or  b) 

2 

II-2 

5-10 

One  plastic  zone, 
free  field  elastic 

4-2 

XX  \€X  Ui 

TTo 

3 

IIa-3 

5-11 

Two  plastic  zones, 
free  field  elastic 

4-3 

lid 

Tib 

3 

IIb-3 

5-12 

Two  plastic  zones, 

(free  field  yielded) 

4-4 

XX  u 

II  (a  orb) 

4 

II-4 

5-13 

Three  plastic  zones, 
(free  field  yielded) 

4-5 
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Table  5-5.  Procedure  for  determining  applicable  solution  branch. 


Case  la 


External  Pressure,  pk 

Phase 

Branch 

Solution  Table 

pb£i[(N  +  l)p,  +  o„] 

1 

1-1 

Table  5-6 

pb>i[(N+l)p.  +  oJ 

2 

1-2 

Table  5-7 

Case  lb 

External  Pressure,  Pv, 

Phase 

Branch 

Solution  Table 

p„<i[(N  +  l)p.+ou] 

1 

1-1 

Table  5-6 

i[(N  +  l)p,+<J„]<p6S 

<*u 

Table  5-7 

2[l-(N+l)v] 

2 

1-2 

of,  Zr  7CT  <  Pb  <  Pb 
2[l-(N+l)vJ 

3 

Ib-3 

Table  5-8 

Pb^Pb 

4 

Ib-4 

Table  5-9 

Case  Ha 

External  Pressure,  ph 

Phase 

Branch 

Solution  Table 

p  <  p*  +  a“ 

Pb  -  2(1- Nv) 

1 

II-l 

Table  5-6 

2(l-Nvj<Pb"Pb 

2 

II-2 

Table  5-10 

P'b  <  Pb  <  Pb 

3 

IIa-3 

Table  5-11 

n.  >p. 

4 

II-4 

Table  5-13 
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Table  5-5.  Procedure  for  determining  applicable  solution  branch  (Continued). 


Case  lib 


External  Pressure,  pb 

_  <  Pa+°u 

Pb  "  2(1- Nv) 

Pa  +^ut  <  Pb  <  Pb 
2(1- Nv) 


Phase 

1 

2 

3 


Branch  Solution  Table 

XI- 1  Table  5-6 

H-2  Table  5-10 

Hb-3  Table  5-12 


Pb^Pb<Pb 

Pb  *  Pb 


4 


II-4 


Table  5-13 


Table  5-6.  Procedure  for  elastic  conditions  (all  cases). 


The  strain  and  stress  fields  are  given  by 

a2 

2Ger  =  -(l-2v)pb  +(pb  -pa)^j 

2Gee  =  -(1-  2v)pb -(pb-Pa)^  =  "2(1- 2v)Pb  -  2Ger 

=-Pb+(Pb-Pa)p- 

a2 

<*0  =  -Pb  “  (Pb  "  Pa  )  =  “2Pb  “  °r 

<*z  =  ~2vPb 
Tunnel  Closure 

^  =  -ee(a)  =  -^[p.-2(l-v)pb] 


(2.11) 

(2.12) 

(2.8) 

(2.9) 

(2.10) 

(2.12’) 
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Table  5-7.  Procedure  for  Case  I  (a  or  b),  Phase  2  solution. 


Compute  R  from 

i 

r  2  (N-l)pb  +  ou~|N~1 

R_a[N  +  l'(N-l)pa+au_ 
Plastic  Zone  Response  (a  <  r  <  R) 

N-l  _ 


(3.28b') 


CL 

(r\ 

=  ~  pn  +  —  -7 

— 

Pa  N-l 

Va  ) 

N-l 


oe  =  Nar-au 
az  =  v(ar+ae) 

2Ge<e)  =  [N  -  (N  +  l)v]ar  -  (1  -  v)au 

2Ge>'  =  -LNmI(n  V^(N ' 1)p- +  a‘ 

2Gee  =  2Ge(ee)  +  2Ge(0p) 

2Ger  =[l-(N+l)v]ar+vau-M2Ge^p) 

Plastic  Zone  Response  (R  <  r  <  °°) 

p'.  =  -°,(R)  =  [p.+^  (f)  "nTI 

r2 

2Ger  =  -(1-  2v)Pb  +  (pb  -  P*a)-r 


N-l 


M+l 


\N-1 


(3.10) 

(3.8) 

(3.19) 

(3.20b) 

(3.24') 


(3.17b) 
(3.15a),  (3.20a) 


2Gee  =  -(1-  2v)pb  -  (pb  -  P*a)^T  =  -2(1-  2v)pb  2Ger 
/  MR2 

Or=-Pb+(Pb-Pa)-2- 


o.—P.-fc-Pw)*' =-2P‘ 

az  = -2vpb 


—  Gr 


(3.10') 

(2.11') 

(2.12') 

(2.8') 

(2.9') 

(2.10’) 
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Table  5-7.  Procedure  for  Case  I  (a  or  b),  Phase  2  solution  (Continued). 

Tunnel  Closure 

^  =  -Ee(a)  =  -^[{N-(N  +  1Hp.-(l-vK] 

1  r(N  +  l)(l-v);(N_i)p  +<T  jlf-T*”  - 1|  (3.20b),  (3.24),  (3.17b) 
2G|_  M  +  N  u  '  JJ|UJ 
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Table  5-8.  Procedure  for  Case  lb,  Phase  3  solution. 


Determine 


boundary  radii  R,  R,  and  R,  and  coefficients  Ai  and  A2. 


Compute  R  from 


(l-2v)cru _ 

R  =  a|  {i_  (N  + 1) v}{(N  -  l)pa  +  <v} . 


N-l 


For  pb  <  pb,  compute 

fB»l  _  1 


B,  I  Crl-Cr2 


1  -C 


2(l-Nv)Apb| 
-1  Crl  J\2(N-l)vApbJ 


and  solve  the  following  nonlinear  equation  for  the  ratio  R . 

(NCrI-Cel)B1(|J,",+(NC,2-Ce2)B2(|J,’",-N(l-2v)pb=0 


(3.33') 


(3.52) 


(3.55) 


Compute  —  from 

a 


_nn-i  N(l-2v)pb-(N-l)CrlB1 


v  a  y 


VRy 


-  (N  -  l)Cr2B2 


^rTY2_1 

R 


(N-l)Pa+au 


(3.58b) 


then  R  =  a 


'Rf1 


N1  R 

and  R  =  -p=p 
K 


R 


Compute  coefficients  from 


(3.52') 

(3.52') 
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Table  5-8.  Procedure  for  Case  lb,  Phase  3  solution  (Continued). 


For  the  case  where  pb  =  pb,  R  °°  an(^ 


_  N(l-2v) 

2  NCr2-C62Pb 


fVT 

<a> 

A1  =  0 


N(l-2v)pb  -(N-  l)Cr2B2 
(N-l)Pa  +  au 


(3.80) 

(3.81b) 


A  ==^- 

2  r_Y2_i 


(3.80’) 


If  pb  <  pb,  the  pressure  at  the  internal  edge  of  the  elastic  zone,  pa  <Jr(R),  is 
required. 

ar(R)  =  CyjAjR71'1  +  Cr2A2R_Y2_1  -  pb  (3.45a’) 

P*a=-tfr(R) 

Inner  plastic  zone  strains  depend  on  circumferential  strains  at  r  -  R  and  r  -  R  so 
those  are  computed  next. 

2Gee  (R)  =  AjR^1"1  +  A2R-Ti_1  +  A  (3.44a’) 

>RfN_l 


/~\  MN  +  2-2(M  +  N+  l)v r/XT  -x  ^  i 
2GME)  ■  (M  +  N)(N  -  l)(l+v)  [  N  ~ 1  P* 


(M  +  2)(l-2v)cu 


(N-1)(M+  l)(l+v) 


-1 


rJ 

+  2Ge0(r) 


^r^n1 


va; 


^  f  r^m+1 


R 


(3.68') 


The  plastic  component  of  £0(r)  is  also  required  in  the  inner  plastic  zone.  It  may  be 
determined  as  follows 

l-2v 


2G£(ep)  (r)  =  2Ge0  (r)  -  2G£(0e)  (r) 


(3.66'),  (3.30') 


(3.17b') 

With  the  information  developed  as  described  above,  stresses  and  strains  at  any 
point  may  be  readily  determined.  For  a  point  at  any  radius  r  >  a,  compare  r  to 
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Table  5-8.  Procedure  for  Case  lb,  Phase  3  solution  (Continued). 


R,  R,  and  R  to  determine  its  location;  then  proceed  to  use  the  appropriate  formulas 
below. 


Tnnpr  Plastic  7,nr»e  Response  (a  <  r  <  R) 


cr=- 


Pa+' 


nN-1  « 

-I 

aj  N- 1 


N-l_ 

ae  =  Nor  -au 
CTz  =  v(or+CJe) 

2Ge(re)  =  [  1  -  (N  + 1)  v]or  +  vou 
2Ge(ee)  =[N-(N+l)v]ar-(l-v)au 

,  ~  \M+1 

2G£?>  =  2Ger(R)[7j 

_(N+1)(1-v)[(N-1)  ]fR 

M  +  N  LV  '  \a 


-  \N-1 

J 


\TJ 


T_ 

R. 


\N-1 


2Gee  =  2Ge(0e)  +  2Ge(0p) 
2Ge(rp)  =  -M2Ge(0p) 
2Ger  =  2Ge(re)  +  2Ge(rp) 


MiHrllp  PI as^>  Zone  Response  (r  <  r  <  R) 

"  s  —  \  M+N  ✓  x  1 

MN  +  2-2(M  +  N  +  l)vrfM_lV  I  -1  -  | 

2G€e  (M  +  N)(N - 1)(1+  v)  "  a  u[lrJ  >a' 


(M  +  2)(l-2v)au 


[rir.il 

+  2Ge.(R)- 

'll' 

y) 

<  F  J 

(N-l)(M  +  l)(l+v) 

MN  +  2  -  2(M  +  N  +  l)v  M-2v 

2Ger  =  -M2Ge0  + - ^ 1+\T  u 


...  "■ 


N-l 


N-l 

-I  + 

a 


N-l 


(3.10) 

(3.8) 

(3.19) 

(3.15a) 

(3.20b) 


(3.70) 

(3.17b) 

(3.15a’) 

(3.17a) 


(3.68) 

(3.69) 
(3.57) 
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Table  5-8.  Procedure  for  Case  lb,  Phase  3  solution  (Continued). 


Og  —  Nhr  Ou 

(3.3a’) 

(3.3a),  (3.3b) 

<?z=CTr 

Outer  Plastic  Zone  Response  (p  <  r  <  R) 

2Ger  =  YiAxr71'1  -  72A2r"Y2_1  +  A 

(3.44a) 

2Gee  =  A^1"1  +  A2r'T2_1  +  A 

(3.44b) 

or  =  C^A^1-1  +  Cr2A2r"Y2_1  -  pb 

(3.45a) 

<re  =  CelAjrYl  +  Ce2A2r  2  -  pb 
oz  =  ^(ae  +  au) 

(3.45b) 

(3.4‘) 

Elastic  Zone  Response  (R  <  r  <  ~) 

R2 

2Ger  =  -(l-2v)Pb  +  (pb  -p*a)^r 

(2. 11') 

2Gee  =  -(1-  2v)pb  -  (pb  -  P*a)^r  =  ~2(1_  2v)Pb "  2Ger 

(2.12') 

,  *1R2 
crr  =  -Pb+(Pb-Pa)^F 

(2.8') 

R2 

Oe  =  -Pb  -  (pb  -  pi  )  =  _2Pb  “  °r 

az  = -2vpb 

(2.9') 

(2.10') 
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Table  5-8.  Procedure  for  Case  lb,  Phase  3  solution  (Continued). 


Tunnel  Closure 


Compute  R,  R,  and  R,  followed  by  Ax  and  A2  as  shown  at  the  beginning  of  the  table. 


Now  compute  circumferential  strain  components  as  follows 
2G€6(r)  =  A1RY‘‘1  +  A2R'Y2"1  +  A 


(3.44a’) 


®M+N 

. .  .  ’T*' 

r  —  m+i  1 

(M  +  2)(l-2v)qu  fRl  _!  +2Ge,(R)-  | 
(n-i)(m+i)(i+v)[1,rJ  J  vR; 

2Ge(0e)(R)  =  ~1_^N+1^va“ 

2Ge(ep)(R)  =  2Ge0  (r)  -  2Ge0e)  (r) 

2Ge0e)  (a)  =  -[N  -  (N  + 1)  v]pa  -  ( 1  -  v)au 

2Gei’> (a)  =  -(N  KN~1)p-  +°"]  _(f 


-nM+1 


(3.68’) 


(3.66’),  (3.30’) 


(3.17b’) 


(3.20b’) 


_  nM+N 


-1  (3.70’) 


Then 


^  =  -e9  (a)  =  -  ^[2Ge?>  (a)  +  2Gei">  («)] 


(3.17b),  (3.66’),  (3.30’),  (3.70’) 
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Table  5-9.  Procedure  for  Case  lb,  Phase  4  solution. 


Determine  boundary  radii  R ,  R ,  and  R ,  and  coefficients  Ai  and  A2. 


Compute  R  from 


R  =  a 


(l-2v)om 


{l-(N  +  l)v}{(N-l)p.  +  a.} 


N-l 


Compute 

(b; 

B, 


crlce2-c„c, 


r2^61  L 


cn,  -c 


'62 

-C 


t2 

61  Crl 


fApb' 

lAPb. 


•  R 

and  solve  the  following  nonlinear  equation  for  the  ratio  ^ . 


(NCrl  — Cei)B1 


^rV1"1 

vRy 


+  (NCr2-C,2)Bl|j  -N(l-2v)pb  =0 


Compute  —  from 


then  R  =  a 


—  n^n-i  N(l-2v)pb  — (N  l^C^Bj  ^ 


Rvrl 


•  (N  -  l)Cr2B5 


^rTY2‘ 

rJ 


(N-l)Pa  +au 


1 

N-l 


R 


and  R  = 

R 

Ir  j 


Compute  coefficients  from 

A  -- 
1  RYl_1 

A,= 


R~Y2_1 


(3.33’) 

(3.84) 

(3.86) 

(3.87b) 


(3.84') 

(3.84') 
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Table  5-9.  Procedure  for  Case  lb,  Phase  4  solution  (Continued). 


Inner  plastic  zone  strains  depend  on  circumferential  strains  at  r  -  K  and  r  -  E, 
while  far  outer  plastic  zone  strains  depend  on  strains  at  r  =  K,  so  those  are 

computed  next. 

2Ge,(R)-AIR"-I  +  A2E-'--‘*A 

2GeW(R)  =  2GE9(R)-2GeW(R) 

2Gee(R)  =  AjR71'1  +  A2R'Yirl  +  A 


(3.44a’) 

(3.88a'),  (3.5a') 

(3.17b') 

(3.44a’) 


MN  +  2-2(M  +  N  +  l)v r/lvT  ^  ^  ] 
2Ge«(R)  =  (M  +  N)(N - !)(l+v)~  N ~ 1  Pa  +  Gu 


(M+N)(N-l)(l+v) 
(M  +  2)(l-2v)an 


-l 


(N  -  1)(M+ 1)(1+  v) 


^R^+1 

VR; 


-1 


vR ) 

2Ge,(R) 


-  \N-1 

ra 

a 


^ipM+1 

vRy 


(3.68') 


The  plastic  component  of  e9(r)  is  also  required  in  the  inner  plastic  zone.  It  may  be 
determined  as  follows 

2G£e1%"l.(N+1)V°u 


(3.66’),  (3.30') 


2G4p)  (r)  =  2Gee  (r)  -  2Ge(ec)  (r) 


(3.17b') 


With  the  information  developed  as  described  above,  stresses  and  strains  at  any 
point  may  be  readily  determined.  For  a  point  at  any  radius  r  >  a,  compare  r  to 
R ,  R,  and  R  to  determine  its  location;  then  proceed  to  use  the  appropriate  formulas 

below. 
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Table  5-9.  Procedure  for  Case  lb,  Phase  4  solution  (Continued). 


Inner  Plastic  Zone  Response 


<*r 


gU 

N-l 


JU  ) 


N-l 


+ 


N-l 


ae  =Ncr-cu 
az  =  v(ar+ae) 

2Ge(re)  =[l-(N  +  l)v]ar  +  vau 
2Ge(0e)  =  [N  -  (N  +  l)v]ar  -  ( 1  -  v)au 


(3.10) 

(3.8) 

(3.19) 

(3.15a) 

(3.20b) 


2Ge(0p)  =  2Ge0p)(R)  - 


_  \M+1 

R 


(N+1)(1_v)j(N-i)  +g  j 

M  +  N  ’  J 


^rnN_1 


(x>\u+l  /.V*-1 


R 


,R 


2Ge0  =  2Ge0e)  +  2Ge0p) 

2Ge(rp)  =  -M2Ge0p) 

2Ger  =  2G£(re)  +  2Ge(rp) 

Middle  Plastic  Zone  Response  (r  <  r  <  R  j 

=  MN  +  2-2(M  +  N  +  l_K  }  ] 

6  (M  +  NKN-l)  1+v  LV 


(M  +  2)(l-2v)cu 


r 


-1 


^  j 

—  xM+1 


N-l 


+  2Ge9(R)(^ 


(N-l)(M  +  l)(l+v) 

MN  +  2-2(M  +  N+l)v  M-2v 

2G£,  =  -M2Ge@  + - ^ 1  T7^Tct“ 


’.-['•**10“** 


Oq  Nar  au 

<Jz=Or 


(3.70) 

(3.17b) 

(3.15a’) 

(3.17a) 


(3.68) 

(3.69) 

(3.57) 

(3.3a') 
(3.3a),  (3.3b) 
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Table  5-9.  Procedure  for  Case  lb,  Phase  4  solution  (Continued). 


Outer  Plast^  Zone  Response  (R<r<R) 

2Ger  =  Y1A1rYl-1-72A2r"Y2"1  +  A 
2Gee  =  A^1-1  +  A2r_Y2_1  +  A 
ar  =  CrlA1rYl_1  +  Cr2A2r  Y2  -  pb 
ae  =  CejA^1"1  +  Ce2A2r"Y2_1  -  pb 

az=^r(oe+cu) 


(3.44a) 

(3.44b) 

(3.45a) 

(3.45b) 

(3.4') 


Far  Outer  Plastic  Zone  Response 
<7r  =  -pb 
<*e  =  “Pb 

c,  =  i(-pb+o„) 

_  [(N+ l)v-N]pb  -vq^ 
zue'  _  (l+v)N 

f(N  +  l)v-Nlpb-vgu 
(l+v)N 


(R  <  r  <  «o) 


2Ge(0e)  = 


rn2 


1- 


2Ge(ep)  =  2Gr£oP)  (r)|^— J  +GC 
2Ge(rp)  =  2GC-2G£(ep)  =  -2Ge^p)(R)f  *  1  +GC 


1+ 


(3.82) 

(3.82) 

(3.5a),  (3.5b) 
(3.88a’),  (3.5a’) 
(3.88b’),  (3.5a’) 

(3.100a’),  (3.100b) 
(3.95),  (3.100a’),  (3.100b) 


2Ge0  =  2Ge0(e)  +  2Gee(p) 
2Ger  =  2Ger(e)  +  2Ge0(p) 


(3.17b’) 

(3.17a’) 
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Table  5-9.  Procedure  for  Case  lb,  Phase  4  solution  (Continued). 


Tunnel  Closure 


Compute  R,  R,  and  R,  followed  by  Ai  and  A2  as  shown  at  the  beginning  of  the  table. 


Now  compute  circumferential  strain  components  as  follows. 

2Gee  (R)  =  A^71'1  +  A2R'72_1  +  A 


(3.44a’) 


vr-MV 

R  ;  i  ^  , 


(M  +  2)(l-2v)an  [RV  _1 
(N-l)(M+l)(l+v)  UJ 


•2Ge,(R){| 


2Gggc)  (r)  = 

2Gegp)  (r)  =  2Ge0  (r)  -  2Ge‘e)  (r) 
2Ge(ee)(a)  =  -[N-(N  +  l)v]pa  -(l-vjo. 


(3.68') 


(3.66'),  (3.30’) 


(3.17b') 

(3.20b') 


Then 


M+l  f  ~ 

2Ger(a)  =  2Ger(R)[f]  -g(^pKN-1>p-'l'q"l  (f)  -1  (3'70,> 


.  =  _e  (a)  =  _  J-[2Ge<;’(a)  +  2Gei'>(a)]  (3.17b),  (3.66'),  (3.30'),  (3.70') 

0 fr  ^ 
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Table  5-10.  Procedure  for  Case  II,  Phase  2  solution. 


For  pb  <  pb  compute 


Determine  plastic  boundary  radius 


the  ratio  (a/R). 


2(l-Nv)APb1  (4.34) 

2(N  -  l)vApb  J 

R  by  solving  the  following  nonlinear  equation  for 


CAI I 


V‘- 


+  Cr2B2|  £ 


a 


-Y2-1 


=  Pb  Pa 


Then  R  =  a/(a/R). 


(4.36) 


Compute  constants 

B,  (4.34') 

Ai  E„-i 

A  =  _A_  (4-34,) 

2  E’*'1 


For  Case  lib  and  pb  equals  pb, 
R  — »  00 

Aj  =  0 


A  —  Pb 

Aa  "  Cr2a"Y2"x 


(4.49b) 


With  the  information  developed  as  described  above,  stresses  and  strains  at  any 
point  may  be  readily  determined.  For  a  point  at  any  radius  r  >  a,  compare  r  to  B  to 
determine  its  location;  then  proceed  to  use  the  appropriate  formulas  below. 
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Table  5-10.  Procedure  for  Case  II,  Phase  2  solution  (Continued). 


Plastic  Zone  Response  (a  <  r  <  R) 


2Ger  =  YiA1riri'1-Y2A2r‘Y2‘1  +  A 

(4.26a) 

2Gee  =  A^r71-1  +  A2r'Y2_1  +  A 

(4.26b) 

cr  =  CrlA1rYl_1  +  Cr2A2r~Y2_1  -  pb 

(4.27a) 

ae  =  CgjAjr"1-1  +  Ce2A2r~Y2_1  —  pb 

(4.27b) 

(4.5b) 

<>z=^(°e  +  °u) 

Elastic  Zone  Response  (R  <  r  <  °°) 


Pi  =  ~or(R)  =  2(1-Nv)pb  +ou 

(4.29) 

R2 

2Ger  =  -(1-  2  v)pb  +  (pb  -  Pi  )^r 

(2.11’) 

2Gee  =  -(1-  2v)pb  -  (pb -P*a)^-  =  -2(1-  2v)pb  -  2Ger 

(2.12') 

(  *1R2 

GT  =-Pb+(Pb-Pa)^2- 

(2.8') 

R2 

<5e  =  -Pb  -  (Pb  -  Pl  )  =  "2Pb  “ 

(2.9') 

cz  = -2vpb 

(2.10’) 

Tunnel  Closure 

^  =  -e!(a)  =  -i(Aia--'  +  A2a--  +  A) 

(4.26b') 
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Table  5-11.  Procedure  for  Case  Ila,  Phase  3  solution. 


Determine  plastic  boundary  radii  K  and  R  and  coefficients  A,  and  A2.  For  pt  <  p„, 
compute 

(4.34) 


[Bxl 

IB, 


crl-cr2 


1  _Cr2l  f  2(l-Nv)Apbj 
-1  Crl  J\2(N-l)vApbJ 


R 

and  then  solve  the  following  nonlinear  equation  for  the  ra  10  R . 

(NC,I-CH)B{|},'I  +  (NC,!-CM)B2(|J,",-N(l-2v)pb=0 


Then  solve  for  the  ratio  R  /  a  from 


N(l-2v)pb-(N-l)CrlBj 


(  rV1'1 

yjt/ 


-(N-l)Cr2B2|  — 


rVY2_1 


(N-l)Pa+cu 


(4.51') 


(4.55b) 


Then 


R  =  a 


rT' 


N-l  R 

and  R  =  7= 


(R/R) 


Compute  constants 

A  =  Bl 

1  RYl_1 

(4.34’) 

A  -  B* 

A2  r-Y2-i 

(4.340 

For  pb  =  pb,  then  R  first  compute 

N(l-2v)pb 

2  NCr2  -  C62 

(4.670 

and  compute  R  from 

f  R Y'1  N(l- 2v)pb  -  (N  -  l)Cr2B2 

UJ  (N-l)p.  +  o. 

(4.680 

A1  —  0  and  A2 
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Table  5-11.  Procedure  for  Case  Ha,  Phase  3  solution  (Continued). 


The  pressure  at  the  internal  edge  of  the  elastic  zone,  pa  is  required  when  pb  <  pb  . 
Therefore  we  compute  the  stress  at  r  =  R  using  the  solution  from  the  outer  plastic 
zone. 


p*a  =  -cr(R)  =  -(CrlA1RYl'1  +  Cr2A2R"Y2'1  -pb) 


(4.27a') 


Inner  Plastic  Zone  Response  (a  <  r  <  R) 

2Gee(R)  =  AjR"1'1  +  A2R'Y2_1  +  A 

MN  +  2-2(M  +  N  +  l)vr,XT  ^  ^  -| 
2GC“=  (M  +  N)(N-l)(l+v-[(N~1)P-+q“1 


M+N 

M 

-1 

UJ 

vN-1 


(M  +  2)(l-2v)qu 


(N-l)(M  +  l)(l+v) 


^R^+1 


\T  J 


-1 


R 


vM+1 


+  2Ge,(E)-  - 


2Ger  =  -M2Ge0  +  1+v 

r  nN-1 


MN  +  2-2(M  +  N  +  l)vg  M-2vg 


1+v 


N-l 


ae  =Nar-au 

<*z=<*r 


(4.26a') 


(4.65') 

(4.66) 

(4.54) 

(4.3a) 
(4.3a),  (4.3b) 


Outer  Plastic  Zone  Response  (R  <  r  <  R) 

2Ger  =  YiA^1-1  -  Y2A2r-Y2-x  +  A 
2Gee=A1rYl-1  +  A2r"Y2_1  +  A 
or  =  CrlA1rYl-1  +  Cr2A2r  Y2  x-pb 
ae  =  CeAr1^1  +  Ce2A2r"Y2_1  -  pb 

<*z=^K+a») 


(4.26a) 

(4.26b) 

(4.27a) 

(4.27b) 

(4.5b) 
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Table  5-11.  Procedure  for  Case  Ha,  Phase  3  solution  (Continued). 


FI  a stic  Zone  Response  (R  <  r  <  ~) 

R2 

2Ger  =-(l-2v)pb  +  (pb-P*a)^r 

(2.11') 

2Gee  =  -(l-2v)pb  -(pb  _p!)“r  =  -2(l-2v)pb  -2Ger 

(2.12') 

(  *)R2 
cyr  =  -Pb+(Pb-Paj^r 

(2.8') 

Oe  =  -Pb  -  (Pb  -  p!  ) “T  =  _2Pb  “ 

(2.9') 

az  =  -2vpb 

(2.10') 

Tunnel  Closure 

Compute  boundary  radii  R  and  R,  and  coefficients  Ai  and  A2  as  shown  at  the 
beginning  of  the  table. 


Then 


and 


2Ge6(R)  -  AjR71"1  +  A2R'T2'1  +  A 


(4.26a) 


AD  MN  +  2-2(M  +  N  +  l)vrfM_1)  +a  1 

D  "  £e^  2G(M+N)(N-1)(1+ v) 


'r^m+n 


-1 


(M  +  2)(l-2v)qu 
+  2G(N-1)(M  +  l)(l+v) 


l)  - 


/-xM+i 

-e,(E)-l®j  (4-65') 
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Table  5-12.  Procedure  for  Case  lib,  Phase  3  solution. 


Determine  plastic  boundary  radius  R,  and  coefficients  Ai  and  A2-  Begin  by 
computing 


JBJ _ 1  f  Ce2  -Cr2ljApb| 

[B2  J  crlce2  -  cr2cei  L-Cei  Crl  J1  APb  J 


and  then  solving  the  following  nonlinear  equation  for  the  ratio 


CmBiIjJ  +Cr2B2 


-Pb  +  Pa  =0 


(4.70) 


(4.72) 


then 

A, 

R1 

a 

{ 

.rJ 

and 

A  — 

(4.70’) 

RTl_1 

A  — 

b2 

(4.70') 

a2  — 

1 

K5 

1 

(-* 

Far  outer  plastic  zone  strains  depend  on  the  circumferential  plastic  strain  at  r  R. 


2Ge,(E). 

=  A1RYl-1  +  A2R-Y2"1  +  A 

(4.26b’) 

2Ge',,(El 

f(N  +  l)v  -  N]pb  -  vou 

(4.73),  (4.16b) 

N(l+v) 

2Ge<,,|(E' 

|  =  2Gee  (e)  -  2Ge(ee)  (e) 

(4.13b) 

With  the  information  developed  as  described  above,  stresses  and  strains  at  any^ 
point  may  be  readily  determined.  For  a  point  at  any  radius  r  >  a,  compare  r  to  R  to 
determine  its  location;  then  proceed  to  use  the  appropriate  formulas  below. 
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Table  5-12.  Procedure  for  Case  lib,  Phase  3  solution  (Continued). 


Tnnpr  Flaatfr  Zone  Response  (a<r<R) 

2Ger  =  YiA/1-1  -  YzAar"72'1  +  A 
2Ge0  =  A/'1-1  +  A2r"T2_1  +  A 
Cr  =  1  +  Cr2A2r  ~  Pb 

ore  =  C^Ajr*1-1  +C62A2r  Y2  1  -pb 


(4.26a) 

(4.26b) 

(4.27a) 

(4.27b) 

(4.5b) 


Far 


Outer  Plastir  Zone  Response  (R  <  r  <  °°) 


Gt  =  “Pb 
<5e  =  "Pb 


°z  =^("Pb+au) 


2Ge(re)  = 


2Ge0e)  = 


[(N  +  l)v-Nlpb-vgu 
(1+ v)N 

[(N  +  l)v-N]pb-vgu 
(1+ v)N 


Rx2 


2Ge<p)  =  2Ge(ep)(R)  -  +GC 


1- 


v  r  y 


2Ge(p)  =  2GC  -  2Ge(p)  =  -2Ge(0p)  (r)I  7  I  +  GC 


1+ 


^2 


2Ge0  =  2Ge0(e)  +  2Ge0(p) 
2Ger  =2Ger(e)  +  2Ge0(p) 


(4.73a) 

(4.73b) 

(4.73c) 
(4.73),  (4.16a) 
(4.73),  (4.16b) 

(4.76'),  (4.77b) 

(4.74),  (4.76'),  (4.77b) 

(4.13b') 

(4.13a’) 
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Table  5-12.  Procedure  for  Case  lib,  Phase  3  solution  (Continued). 


Tunnel  Closure 


Compute  R ,  Ai,  and  A2  as  shown  at  the  beginning  of  the  table. 

TT  =  '£e(a)  =  + + A1 


(4.26b') 
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Table  5-13.  Procedure  for  Case  II,  Phase  4  solution. 


Determine  plastic  boundary  radii  R  and  R.  Begin  by  computing 


'*1]  _  1 

iBj  CrlCe2 -Cr2Cei 


'02 


-C 


r2 


-C91  Crl  j 


'Apb 

LAPb. 


R 


and  then  solving  the  following  nonlinear  equation  for  the  ratio  ^ . 


(NCrt-C.jBi 


vRy 


•  (NCr2  -  Ce2  )B2 


^rVY2_1 

vRy 


-N(l-2v)pb  =  0 


Then  solve  for  the  ratio  R  /  a  from 

_  N_!  N(l- 2 v)pb  -  (N -  l)CrlBx 


Then 


f  R 


R  =  a 


'lpYl_1 


-(N-l)Cr2B2 


^Rv7rl 

A 

R 


(N-l)Pa+au 


r^n  1 


n-i  -  R 
and  R  =  /_ 


(e/k) 


Compute  constants 
Al  RYl_1 


A9  — 


b2 


R 


-Y2-1 


(4.70) 


(4.87) 


(4.88b) 


(4.70) 

(4.70) 


Far  outer  plastic  zone  strains  depend  on  the  circumferential  plastic  strain  at  r  -  R. 
2Gee  (R)  =  A.R^1  +  A^2'1  +  A  <4-26b,) 

on(e)/p\  [(N+l)v-N]pb-yqa_  (4.73),  (4.16b) 

2Ge®  lRr  n(i+V) 

2Ge(ep)  (r)  =  2Gee(R)-2Ge(ee)(R)  (4-13b) 

With  the  information  developed  as  described  above,  stresses  and  strains  at  any^ 
point  may  be  readily  determined.  For  a  point  at  any  radius  i  r  >  a  compare  r  to  R  to 
determine  its  location;  then  proceed  to  use  the  appropriate  formulas  below. 
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Table  5-13.  Procedure  for  Case  II,  Phase  4  solution  (Continued). 


Trmor  Plastic  Zone  Response  (a<r<R) 

2Ge0  (R)  =  AaRYl-1  +  A2R'Y2_1  +  A 

MN  +  2-2(M  +  N+1)v[(TvT_1)  i 
2G£e_  (M  +  N)(N-l)(l+v)  11  jPa  uJ 


R 

r 


M+N 


(M  +  2)(l-2v)au 


(N-l)(M  +  l)(l+v) 


v  M+l 


- 1  -1 


+  2Ge,(R)|f 


-l 

N  M+l 

J 


N-l 


2Ger  =  -M2Ge0  + 

a9  =  Nor-ou 
Oz=Gr 


MN  +  2  -  2(M  +  N  +  l)v  M-2v 


1+v 


Gr  .  vu 
1+V 


\N-1 


( rY  '  .  G„ 


+ 


N-l 


(4.26a’) 


(4.65’) 

(4.66) 

(4.54) 

(4.3a) 
(4.3a),  (4.3b) 


Outer  Plastic  Zone  Response  ^R  <  r  <  R 

2Ger  =  YxA^1"1  -  Y2A2r_Y2"1  +  A 
2G£e  =  Ajr^  +  Aar-^  +  A 
Gr  =  CrlAxrYl  1  +  Cr2A2r  72  1_  Pb 
ce  =  CelA1rYl-1  +  Ce2A2r  72  1_pb 

az=^r(oe  +  ^u) 


(4.26a) 

(4.26b) 

(4.27a) 

(4.27b) 

(4.5b) 
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Table  5-13.  Procedure  for  Case  II,  Phase  4  solution  (Continued). 


Far  Outer  Plastic  Zone  Response 


(R  <  r  <  °°) 


ar  =  -pb 
Gq  =  —  Pb 

o,  =  i(-pk  +  o„) 

-c  ,.,  KN+i)v-N]pt-yq. 

20e'  "  (l+v)N 

_[(N  +  l)v-N]pb-vau 
(l+v)N 


2Geie)  = 


2Ger  =  2Ger(R)(^ 


+GC 


1- 


vr  y 


2Ge(rp)  =  2GC  -  2Ge(ep)  =  -2Ge|,p)(R) 


-I  +GC 


1+ 


2Ge0  =  2Ge0(e)  +2Ge0(p) 

2Ger  =2Ger(e)+2Ge0(p) 


(4.73a) 

(4.73b) 

(4.73c) 
(4.73),  (4.16a) 
(4.73),  (4.16b) 

(4.76'),  (4.77b) 

(4.74),  (4.76’),  (4.77b) 

(4.13b') 

(4.13a') 


Timnel  Closure 

Compute  boundary  radii  R  and  R,  and  coefficients  Ai  and  A2  as  shown  at  the 
beginning  of  the  table.  Then 

2Ge0  (R)  =  A,Rr -1  +  A2R“Y2_1  +  A 

and 

r/5?\M+N 

_ _  ~  A/«  r  TLT  .  1  .  I  (  T?  \ 

-1 


(4.26a) 


AD_  ,  x  MN  +  2-2(M  +  N  +  l)v_[/M_1y  ]| 

“5““  £e(a)  2G(M  +  N)(N  - 1)(1+  v) 


(M  +  2)(l-2v)gu 
+  2G(N-1)(M  +  1)(1+v) 


(Sf1- 


vay 


-£« 


(4.65') 
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SECTION  6 

NUMERICAL  EXAMPLES 

Numerical  comparisons  of  the  results  from  finite  element  calculations  with  the 
results  from  the  closed-form  solutions  discussed  in  the  preceding  sections  are 
presented  in  this  section.  Two  examples  for  each  of  the  four  cases  (Cases  la,  lb,  Ha, 
and  lib)  are  presented:  one  example  with  M  equal  to  N  (associative  flow  rule)  and 
one  example  with  M  equal  to  one  (nonassociative  flow  rule).  Results  using  the  two 
solution  techniques  are  presented  in  terms  of:  (1)  a  table  of  differences  in  computed 
closure  for  a  given  load  level;  (2)  plots  of  external  pressure,  pb,  versus  closure;  and 

(3)  plots  of  stress  and  displacement  versus  normalized  radius  (r/a)  at  given  load 
levels.  The  material  properties  and  loads  used  for  each  example  are  provided  in 
Table  6-1.  The  material  properties  were  chosen  to  correspond  to  a  soft  rock.  For 
each  example,  the  radius  of  the  hole  was  taken  as  unity  and  loading  was  applied 
until  closure  of  approximately  five  percent  was  achieved.  Except  for  the  Case  la 
examples,  where  free-field  yielding  cannot  occur,  the  load  levels  were  high  enough  to 

cause  the  free  field  to  yield. 

Prior  to  performing  the  nonlinear  finite  element  calculations,  a  mesh  discretization 
study  was  performed  using  a  linear  material.  The  objective  of  that  study  was  to 
determine:  (1)  the  size  of  the  mesh  required  to  approximate  an  infinite  medium;  and 
(2)  the  number  of  elements  required  to  achieve  accurate  results.  For  the  linear 
problem,  we  used  Youngs  modulus  of  one  million  psi,  Poisson's  ratio  of  0.25,  and  a 
hole  radius  of  one  inch.  Loading  was  applied  in  two  steps.  In  the  first  step,  the 
internal  and  external  pressures  were  both  equal  to  one  thousand  psi.  In  the  second 
step,  the  internal  pressure  was  held  at  one  thousand  psi,  while  the  external  pressure 
was  increased  to  two  thousand  psi.  A  mesh  containing  four  hundred  elements  with 
an  external  radius  of  two  hundred  inches  gave  displacement  and  stress  results 

accurate  to  four  decimal  places.  This  level  of  mesh  refinement  may  be  overly 

conservative,  but  was  selected  to  eliminate  the  possibility  of  a  coarse  mesh  being  a 
significant  source  of  error  in  the  nonlinear  finite  element  calculations.  The  linear 
analysis  results  are  presented  in  Figure  6-1.  Note  that  the  finite  element  results 
(dashed  lines)  are  not  distinguishable  from  the  closed-form  solution  (solid  lines).  In 
Figure  6-1  and  all  subsequent  figures,  compression  is  shown  positive. 
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Table  6-1.  Material  and  load  parameters  for  example  problems 
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Displacement  (in)  Stress  (psi) 


1200 


800  H 


400  A 


0 


0 


-  Analytical 

_ _ _  Finite  Element 

- - — 1 - 

5 

- ! - 

10 

15  2< 

Normalized  Radius  (r/a) 
a.  Step  1  stress  profiles. 


b.  Step  1  displacement  profile. 

Figure  6-1.  Comparisons  of  analytical  and  finite  element  analyses  for  linear 
material  response. 
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Displacement  (in) 


3200 


2400- 


m  1600- 

I 

m 


—  Analytical 
_  Finite  Element 


5  10 

Normalized  Radius  (r/a) 

c.  Step  2  stress  profiles. 


d.  Step  2  displacement  profile. 

Figure  6-1.  Comparisons  of  analytical  and  finite  element  analyses  for  Unear 
material  response  (Continued). 
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Closure  results  at  the  final  load  level  for  each  of  the  eight  nonlinear  cases  are  pre¬ 
sented  in  Table  6-2.  The  maximum  difference  in  closure  between  the  finite  element 
and  the  closed-form  solutions  was  0.337  percent.  In  performing  these  analyses,  we 
found  that  the  differences  in  computed  closures  could  be  lessened  by  reducing  the 
load  increment  and  tightening  the  convergence  tolerances  in  the  nonlinear  finite 
element  analyses.  This  was  particularly  true  in  the  nonassociative  flow  examples, 
where  the  elastic-plastic  constitutive  matrix  (required  for  implicit  finite  element 
analysis  utilizing  Newton-Raphson  techniques  for  solving  nonlinear  equations) 
becomes  nonsymmetric.  Our  finite  element  program  utilized  a  symmetric  equation 
solver,  and  thus  an  approximation  to  the  true  elastic-plastic  constitutive  matrix  had 
to  be  used.  Although  we  may  have  been  able  to  achieve  even  closer  matches  between 
analytical  and  numerical  results,  the  differences  presented  in  Table  6-2  were  judged 
to  be  small  enough  to  convince  us  that  the  hypotheses  used  in  developing  the 
analytical  solutions  were  justified. 

Plots  of  external  pressure,  pb,  versus  closure  for  the  eight  example  problems  are 
shown  in  Figures  6-2  through  6-9.  There  is  no  discernible  difference  between  the 
closed-form  and  finite  element  results. 

Plots  of  stress  (ar,  ae,  and  oz)  versus  normalized  radius  (r/a)  and  radial  dis¬ 
placement  versus  normalized  radius  (r/a)  for  each  of  the  eight  cases  are  presented  in 
Figures  6-10  through  6-17.  Again,  there  is  no  discernible  difference  between  the 
analytical  and  finite  element  results. 


106 


Table  6-2.  Closure  results. 


Example 

Number 


Applied  External 
Pressure,  pb  (psi) 

Closure  {%) 
(Finite  Element) 

Closure  (%) 
(Analytical) 

Difference  (%) 
* 

1,100 

5.7390 

5.7481 

-0.158 

675 

5.0254 

5.0347 

-0.185 

750 

5.0573 

5.0744 

-0.337 

1,200 

4.9428 

4.9479 

-0.103 

1,500 

5.0083 

5.0222 

-0.277 

930 

4.9811 

4.9976 

-0.330 

1,000 

4.9895 

5.0021 

-0.252 

1,425 

5.1112 

5.1180 

-0.117 

8 


*  Difference  (in  percent)  is  computed  by  subtracting  the  finite  element  solution  from 
the  analytical  ^solution,  then  dividing  by  the  analytical  solution,  and  mulfrplymg 

by  100. 
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Figure  6-2.  Pressure  versus  closure  (Example  1). 


Figure  6-3.  Pressure  versus  closure  (Example  2). 
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Pressure  (psi)  Pressure  (psi) 
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Pressure  (psi)  Pressure  (psi) 


Figure  6-6.  Pressure  versus  closure  (Example  5). 


Figure  6-7.  Pressure  versus  closure  (Example  6). 
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Ill 


Displacement  (in) 


a.  Stress  profile. 


b.  Displacement  profile. 

Figure  6-10.  Example  1  response  profiles. 
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Displacement  (in) 


a.  Stress  profile. 
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Displacement  (in) 


a.  Stress  profile. 


b.  Displacement  profile. 


Figure  6-12.  Example  3  response  profiles. 
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Displacement  (in) 
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Displacement  (in) 


3000 


Analytical 
Finite  Element 


a.  Stress  profile. 


b.  Displacement  profile. 

Figure  6-14.  Example  5  response  profiles. 
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Displacement  (in) 


Displacement  (in) 


b.  Displacement  profile. 

Figure  6-16.  Example  7  response  profiles. 
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Displacement  (in) 


2400 


Analytical 
Finite  Element 


Figure  6-17.  Example  8  response  profiles. 


L9 


SECTION  7 
REFERENCES 


Detoumay,  E.,  C.  M.  St.  John,  and  D.  E.  Van  Dillen,  "An  Investigation  of  the  Failure 
Resistance  of  Rockbolted  Tunnels  for  Deep  Missile  Basing,"  DNA-TR-82-183, 
Agbabian  Associates,  2  September  1983.  (UNCLASSIFIED) 

Florence,  A.  L.,  Jr.  and  L.  E.  Schwer,  "Axisymmetric  Compression  of  a  Mohr- 
Coulomb  Medium  Around  a  Circular  Hole,"  International  Journal  for  Numerical  and 
Analytical  Methods  in  Geomechanics,  Vol.  2,  pp.  367-379,  SRI  International,  1978. 
(UNCLASSIFIED) 

Florence,  A.  L.,  Jr.  and  L.  E.  Schwer,  "Axisymmetric  Compression  of  a  Mohr- 
Coulomb  Medium  Around  a  Circular  Hole,"  SRI  International,  UNPUBLISHED, 
Undated.  (UNCLASSIFIED) 

Hendron,  A.  J.,  Jr.  and  A.  K.  Aiyer,  "Stresses  and  Strains  Around  a  Cylindrical 
Tunnel  in  an  Elasto-Plastic  Material  with  Dilatancy,"  Technical  Report  No.  10,  U.S. 
Army  Corps  of  Engineers,  Omaha  District,  September  1972.  (UNCLASSIFIED) 

Merkle,  D.  H.,  Notes  on  modifying  the  Hendron  and  Aiyer  Solution,  Applied 
Research  Associates,  UNPUBLISHED,  3  September  1982.  (UNCLASSIFIED) 

Newmark,  N.  M.,  "Design  of  Rock  Silo  and  Rock  Cavity  Linings,"  Appendix  II  of 
Ground  Motion  Technology  Review.  SAMSO  TR-70-114,  Nathan  M.  Newmark 
Consulting  Engineering  Services,  August  1969.  (UNCLASSIFIED) 

Wintergerst,  G.  L.,  D.  N.  Burgess,  J.  L.  Merritt,  and  P.  E.  Senseny,  "Axisymmetric 
Compression  of  a  Mohr-Coulomb  Medium  with  Arbitrary  Dilatancy,"  DNA-TR-91- 
38,  BDM  Engineering  Services  Company,  1  November  1991.  (UNCLASSIFIED) 


DISTRIBUTION  LIST 
DSWATR-97-5 


DEPARTMENT  OF  DEFENSE 

DEFENSE  INTELLIGENCE  AGENCY 
ATTN:  TWJ 

DEFENSE  SPECIAL  WEAPONS  AGENCY 
ATTN:  OPS 

ATTN:  PMCO,  MR  D  BRUDER 
ATTN:  PMT,  K  GOERING 
ATTN:  PMT,  P  SENSENY 

2CY  ATTN:  TRC 

ATTN:  WEL,  D  RICE 
ATTN:  WEP,  W  ULLRICH 

DEFENSE  TECHNICAL  INFORMATION  CENTER 

2  CY  ATTN:  DTIC/OCP 

FC  DEFENSE  SPECIAL  WEAPONS  AGENCY 
ATTN: FCTO 
ATTN:  FCTOS 
ATTN:  FCTT  DR  BALADI 
ATTN:  FCTTS,  E  TREMBA 

DEPARTMENT  OF  THE  ARMY 

ARMY  RESEARCH  PROJECT  AGENCY 

ATTN:  DEFENSE  SCIENCES  OFFICE 

U  S  ARMY  CORPS  OF  ENGINEERS 
ATTN:  CERD-L 

U  S  ARMY  ENGINEER  DIST  OMAHA 
ATTN:  MROED-S,  H  GAUBE 

US  ARMY  ENGR  WATERWAYS  EXPER  STATION 
ATTN:  C  D  NORMAN 
ATTN:  CEWES-SD 
ATTN:  CEWES-SE,  MR  L  K  DAVIS 
ATTN:  CEWES-GS,  DR  D  BANKS 
ATTN:  WESGR-M,  J  WARRINER 
ATTN:  TECHNICAL  LIBRARY 
ATTN:  W  MCMAHON 
ATTN:  W  MILLER 

DEPARTMENT  OF  THE  AIR  FORCE 

AIR  FORCE  CTR  FOR  STUDIES  &  ANALYSIS 
ATTN:  AFSAA/SAI,  RM  1D363 
THE  PENTAGON 


AIR  UNIVERSITY  LIBRARY 
ATTN:  AUL - LSE 

DEPARTMENT  OF  ENERGY 

LOS  ALAMOS  NATIONAL  LABORATORY 
ATTN:  REPORT  LIBRARY 

SANDIA  NATIONAL  LABORATORIES 
ATTN:  A  F  FOSSUM 
ATTN:  J  D  ROGERS 
ATTN:  R  ZIMMERMAN 
ATTN:  T  BERGSTRESSER 
ATTN:  TECH  LIB 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

ANALYTIC  SERVICES,  INC.  (ANSER) 

ATTN:  K  BAKER 

APPLIED  RESEARCH  ASSOCIATES,  INC. 

ATTN:  C  J  HIGGINS 

APPLIED  RESEARCH  ASSOCIATES,  INC. 

ATTN:  S  BLOUIN 

APPLIED  RESEARCH  ASSOCIATES,  INC. 

ATTN:  R  FRANK 

APTEK,  INC. 

ATTN:  B  LEWIS 
ATTN:  Y  MURRAY 

BDM  INTERNATIONAL  INC 
2CY  ATTN:  DANNY  N  BURGESS 
2  CY  ATTN:  DAVID  M  KENDALL 
2  CY  ATTN:  GARY  L  WINTERGERST 

BOEING  TECHNICAL  &  MANAGEMENT  SYCS,  INC. 
ATTN:  W  M  LEAVENS 
ATTN:  AW  SPENCER 
ATTN:  R  BRYAN  CAIRNS 

ITT  RESEARCH  INSTITUTE 

ATTN:  DOCUMENT  LIBRARY 
ATTN:  M  JOHNSON 

JAYCOR 

ATTN:  CYRUS  P  KNOWLES 


Dist-1 


DSWA-TR-97-5  (DL  CONTINUED) 


KAMAN  SCIENCES  CORP. 

ATTN:  RICHARD  KEEFFE 

KAMAN  SCIENCES  CORPORATION 
ATTN:  DASIAC 
ATTN:  DASIAC/DARE 

LACHEL  AND  ASSOCIATES,  INC 
ATTN:  C  LINAMEN 
ATTN:  J  BECK 

LOGICON  R  AND  D  ASSOCIATES 
ATTN:  DR  T  A  PUCIK 
ATTN:  LIBRARY 

LOGICON  R  AND  D  ASSOCIATES 
ATTN:  E  HUMPHREYS 

MAXWELL  TECHNOLOGIES,  INC 
ATTN:  K  D  PYATT,  JR 

NTS  ENGINEERING 

ATTN:  S  SHORT 

PACIFIC-SIERRA  RESEARCH  CORP. 

ATTN:  H  BRODE 

PACIFIC-SIERRA  RESEARCH  CORP. 

ATTN:  D  GORMLEY 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  H  PRATT 

ATTN:  TECHNICAL  REPORT  SYSTEM 


SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  W  LAYSON 

SOUTHWEST  RESEARCH  INSTITUTE 
ATTN:  B  THACKER 

SRI  INTERNATIONAL 

ATTN:  DR  JIM  GRAN 

THE  AEROSPACE  CORP 

ATTN.  LIBRARY  ACQUISITION  Ml/199 

TITAN  CORPORATION 

TITAN  RESEARCH  &  TECHNOLOGY  DIVISION 
ATTN:  J  THOMSEN 

TITAN  CORPORATION  (THE) 

TITAN  RESEARCH  &  TECHNOLOGY  DIVSN 
ATTN:  LIBRARY 
ATTN:  S  SCHUSTER 
ATTN:  Y  M ITO 

TRW  SPACE  &  DEFENSE  SECTOR  SPACE  & 
TECH  GROUP 

ATTN:  OUT6,  W  WAMPLER 

UTD,  INC. 

ATTN:  E  FOSTER 

WEIDLINGER  ASSOC,  INC. 

ATTN:  E  WONG 
ATTN:  H  LEVINE 

WEILDLINGER  ASSOCIATES,  INC. 

ATTN:  I  SANDLER 


Dist-2 


